UW/PT 04-15 



Staggered Chiral Perturbation Theory at Next-to-Leading Order 

Stephen R. SharpcQ and Ruth S. Van de WateiQ 

Physics Department, University of Washington, Seattle, WA 98195-1560 
(Dated: February 1, 2008) 

We study taste and Euclidean rotational symmetry violation for staggered fermions at nonzero lat- 
tice spacing using staggered chiral perturbation theory. We extend the staggered chiral Lagrangian 
to O(aV), 0(a 4 ) and 0{a 2 m), the orders necessary for a full next-to-leading order calculation of 
pseudo-Goldstone boson masses and decay constants including analytic terms. We then calculate 
a number of 5*0(4) taste-breaking quantities, which involve only a small subset of these NLO op- 
erators. We predict relationships between 5*0(4) taste-breaking splittings in masses, pseudoscalar 
decay constants, and dispersion relations. We also find predictions for a few quantities that are not 
50(4) breaking. All these results hold also for theories in which the fourth-root of the fermionic 
determinant is taken to reduce the number of quark tastes; testing them will therefore provide 
] evidence for or against the validity of this trick. 
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Simulations of lattice QCD using staggered fermions Q are able, at present, to reach significantly smaller dynamical 
quark masses than those using other fermion discretizations 0, 0. This computational advantage is due in part 
to the axial U(l) symmetry remaining at nonzero lattice spacing, which protects the quark mass from additive 
renormalization. The primary disadvantage of staggered fermions is that each continuum flavor comes in four tastes. 
In the continuum limit, the tastes become degenerate and one can formally remove the extra tastes by taking the 
fourth-root of the quark determinant (the "v / Det trick"). At non-zero lattice spacing, however, the SU(4) taste 
symmetry is broken by discretization errors of C(a 2 ), where a is the lattice spacing. These errors are numerically 
significant in present simulations This has two related implications. First, since one must, in practice, take 
\ the fourth-root of the determinant before taking the continuum limit, the resulting underlying quark action is likely 
non-local. It is thus not guaranteed to lie in the same universality class as QCD, even if the non-local part is, in some 
sense, of 0(a 2 ) and thus small. 1 Second, assuming that the continuum limit is correct, it is necessary to perform a 
combined continuum and chiral extrapolation taking into account taste-breaking effects of 0(a 2 ) 0- 

The first issue is clearly more fundamental, but also more difficult to address. To date, support for the use of 
the \ / T5et trick comes from an accumulation of numerical or indirect evidence: the accurate agreement of lattice 
J^" 1, and experimental results for all "gold-plated" quantities that have been calculated Q ; the success of fits of the light 
pseudoscalar meson masses and decay constants to the predictions of chiral perturbation theory including the 0{a 2 ) 
taste- breaking errors ; the near four-fold degeneracy of eigenvalues of the staggered Dirac operator and its ability 
to correctly measure the topology of the gauge configurations 0, ; and the results of analytic studies of the time 
continuum limit [Tfl lllj . None of these tests purports to be conclusive. Nevertheless, in the absence of a more 
theoretical approach, it is important to have further numerical tests. This is one of the motivations for the present 
?3 i work. 

Our major focus, however, is on the further development of chiral perturbation theory for staggered fermions 
including discretization errors. This "staggered chiral perturbation theory" (or S^PT for short) is the theoretical tool 
which determines the functional forms needed to do the combined continuum and chiral extrapolations incorporating 
taste violations. The development of S%PT for staggered fermions without the v^Det trick, i.e. with four tastes per 
flavor, is a standard application of the methods of effective field theory p^,IT^|. In order to apply SvPT to simulations 
done with the v^Det trick, however, one must alter the X-PT Feynman rules by hand (or alternatively use a variant 
of the rep lica trick .14]), so there is no underlying "staggered chiral Lagrangian" which can reproduce the Feynman 
rules |13| . Nevertheless, if the resulting expressions describe the outcomes of numerical simulations done using the 
V^Det trick, one gains confidence that the long distance physics of the staggered simulations reproduces that of QCD 
in the continuum limit. Indeed, the expressions in S^PT go over to those of standard QCD xPT in the continuum 
limit, and give an explicit example of how the effect of the v'Det trick can vanish smoothly. 
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1 For summaries of this and related problems with the v'Det trick see Refs. 00 
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We consider here only SxPT applied to the properties of pseudo-Goldstone bosons (PGBs). 2 SxPT is a combined 
expansion in powers of the quark masses and lattice spacing, in which the usual power counting is m ~ a 2 . The leading 
order (LO) correction to continuum x?T is the 0(a 2 ) potential, which was determined for a single flavor of staggered 
fermion in Ref. [T^]. This potential retains an 50(4) subgroup of the continuum SU(4) taste symmetry. Since this 
subgroup is larger than the discrete lattice symmetry group, the masses of the PGBs in certain lattice irreps are 
predicted, at LO, to be degenerate. Reference generalized the potential to multiple staggered flavors, and showed 
that the degeneracies remain. These predictions of SxPT work very well; both in old quenched results |15l Ha . Il7j and 
in recent dynamical simulations (using the v'Det trick) jj]. This success gives one confidence in the utility of SxPT. 

It is well known from continuum xPT that practical applications to the properties of the physical PGB mesons 
requires one to work at least to next-to-leading order (NLO). The same is true of chiral extrapolations of lattice 
data |18| . Thus for SxPT to be a practical tool for guiding extrapolations, it must be extended beyond LO. Important 
first steps were taken in Refs. 0] and where the one-loop contributions to the mass and decay constant of the 
lattice Goldstone boson (that with taste £5) were calculated. These forms, and in particular the taste-breaking built 
into them, have proved essential to describe the numerical results 0. Conversely, the success of the fits provides more 
evidence for the applicability of SxPT, and the lack of problems introduced by the \/15et trick. 

In this paper we continue the extension of SxPT to NLO by determining all the operators proportional to a 2 p 2 , 
a 2 m and a 4 , including source terms for left- and right-handed currents and for scalar and pseudoscalar densities. 
These operators give the analytic NLO contributions that incorporate the effects of discretization. Each is multiplied 
by a different, unknown low energy coefficient (LEC). The full NLO result is obtained by adding their contributions 
to those from one-loop diagrams and from the continuum NLO operators. 

In continuum xPT there are of order ten NLO operators, and it is possible to consider enough physical quantities 
to determine their coefficients and make predictions. In SxPT we find that there are an order of magnitude more 
NLO operators. While, in principle, there is a corresponding increase in the number of physical quantities that can 
be calculated (e.g. scattering amplitudes with many choices of external tastes), in practice it will be difficult to carry 
through a NLO program along the lines of that in the continuum. 

Nevertheless, our enumeration of operators is not a purely academic exercise, for there is one set of quantities for 
which very few operators contribute, and a number of new predictions can be made. These are quantities which 
are non- vanishing only because the 50(4) taste symmetry is broken. Since this symmetry is preserved by the LO 
Lagrangian, it is also preserved by the one-loop contributions, as these involve LO vertices. Thus the only source of 
50(4) breaking is the NLO operators. Furthermore, we find that for pion masses and pseudoscalar decay constants 
(i.e. the vacuum to pion matrix element of the pseudoscalar density), the only contribution to 50(4) breaking comes 
from the operators proportional to a 2 p 2 , and that there are few enough of these that we can make several predictions. 
These are valid only at NLO, being violated by higher order terms. Comparison of these SxPT predictions with 
lattice data numerically tests the yDet trick. 

Discussions of staggered fermions can quickly become rather technical, so in an attempt to keep the main text of 
this paper accessible, we relegate most of the details of operator enumeration to an appendix. The precise organization 
is as follows. In Sections HT1 and llTfl we review the construction of the continuum quark-level effective staggered action 
and leading order staggered chiral perturbation theory, respectively. In Section HVI we discuss the symmetry-breaking 
pattern exhibited by the (D(a 2 p 2 ), 0(a 4 ), and 0(a 2 m) operators, and outline the general consequences for various 
PGB properties. Our specific predictions for relations between 50(4) breaking in masses, dispersion relations and 
pseudoscalar decay constants are presented and explained in Section We conclude in Section IVll We determine 
the NLO staggered operators in Appendix 1X1 and collect them in the tables of Appendix IB"1 

II. CONTINUUM EFFECTIVE LAGRANGIAN FOR STAGGERED QUARKS 

The construction of the chiral Lagrangian including discretization errors proceeds in two steps |2(j. One first 
determines the continuum effective Lagrangian for quarks with momenta much smaller than the lattice cut-off 1/a. 
The dominant power-law dependence on the lattice spacing is thereby made explicit. In the second step one maps the 
continuum quark-level Lagrangian onto an effective chiral Lagrangian. In this section we briefly review the first step 
for staggered fermions, including in the effective Lagrangian the leading discretization effects, which are proportional 



2 There is a potential ambiguity in this terminology. We use the term PGB to refer to all pseudoscalar mesons (and fermions, once we 
consider partially quenched theories) whose masses vanish in the combined chiral and continuum limits. This includes the pion whose 
mass vanishes in the chiral limit even at non-zero lattice spacing, which we refer to as the "lattice Goldstone" pion, as well as all the 
other light pseudoscalars with different tastes. 
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to a 2 . This Lagrangian was determined for a single flavor in Ref. [I2(, and generalized to multiple flavors in Ref. |13| . 
While we obtain no new results here, we attempt to clarify aspects of the discussion in these papers. 

The construction of the effective quark Lagrangian follows the standard method developed by Symanzik [^3, m 
which one determines all local continuum operators of dimension up to and including six that are invariant under the 
lattice symmetry group. Given the complexity of the symmetry group for staggered fermions [22^ . a simple way to 
perform this enumeration is to first write down allowed lattice operators of dimension up to six, and then match them 
onto continuum operators. The first step utilizes all of the lattice symmetries, including translations. It is discussed 
in detail in Appendix A of Ref. 0] ■ The second step uses only transformations that map the 2 4 hypercube onto itself, 
which are related to spin and taste transformations in the continuum p3j . As is well known, there are no dimension 
5 operators |24ll25|. At dimension 6 one finds purel y gl uonic operators, fermion bilinears and four-fermion operators. 
These operators were enumerated in Refs. [2(j and [l2j|. 

Neither the gluonic operators nor the fermion bilinears break the continuum taste symmetry because they do not 
contain any taste matrices. Since our major focus is taste symmetry breaking, we do not reproduce these operators 
here. They do, however, break the continuum 5*0(4) rotation symmetry down to its hypercubic subgroup, SW4, and 
we discuss this issue briefly in Appendix El 

SU(4) taste symmetry is broken by the four-fermion operators. We first recall the results for a single flavor and 
then gen eralize to multiple flavors. We write quark bilinears using the direct-product (spin <g) taste) notation of 
Refs. [23, [23, [23 3 Each four fermion operator in the effective Lagrangian turns out to be a product of two bilinears 
with the same spin and taste. Thus one can use a compact notation to represent them. Following Refs. 0, [2^| we 
write the pair of spin matrices in the two bilinears as 

S = 1® 1, P = 7 5®75, ^ = ^7p®7„> 4 = $^7m6<8>76/*, t = ® > U) 

using the definitions 7 M „ = 7^7,, and 7^,5 = 7^75. The pair of taste matrices are denoted in the same manner with 
7^ — > Cm- The four fermion operator is then denoted by the spin and taste matrices that it contains. For example 

[A x T] = ^ 5 ® &p)QQ(^ ® Zp»)Q > ( 2 ) 

where Q, Q are Dirac fields with spinor and taste indices. We stress that Q, Q are at this stage continuum fields, and 
that the two Q's and two Q's in this operator reside at the same position in Euclidean space. 

The fields Q, Q also have color indices, which can be contracted in @ either within each bilinear or between them. 
These two color contractions are related by Fierz transformations. In Ref. jl2j |. these transformations were used to 
select a basis in which all operators had the same type of color contraction. Here we proceed differently, following 
Ref. 0]. We keep both color contractions, and use the notation such as [A x T] to implicitly refer to both. This is a 
sensible notation because, when we map these operators into the chiral Lagrangian, the choice of color contraction is 
irrelevant: the two operators map onto the same chiral operators. Using this notation allows us to shorten the list of 
operators, since by Fierz transformations we can always bring the operators into a form in which both bilinears have 
"odd distance" (see Ref. 0] for an explanation of this terminology). For example, the operator {[S x S] — [P x P]}, 
which consists of two bilinears both of distance zero, can be transformed into a linear combination of operators with 
odd-distance bilinears. 

We now list the resulting four-fermion operators. They can be divided into two types based on how badly they 
break taste symmetry. Twelve have their spin and taste indices contracted separately, as in the example of (|2"|. and 
are therefore rotationally- invariant: 

£ FF(A) ~ [S x A] + [S x V] + [A x S] + [V x S] + [P x V] + [P x A] + [V x P] + [A x P] + 

[T x V] + [T x A] + [V x T] + [A x T] . (3) 

Each operator has a unique coefficient, which we do not show, that is proportional to a 2 , but also depends on g(a) 2 
and ln(a). Because the indices on the taste matrices are contracted between bilinears, these operators preserve an 
50(4) subgroup of the full SU(4) taste symmetry. The four remaining four fermion operators have spin and taste 
indices contracted together: 

< F(B) ~ x V,] + [T, x Afj] + [V, x ly + [A„ x T M ] . (4) 



In the standard staggered convention, the spin matrices, 73, are Euclidean gamma matrices, while the taste matrices, fj*, are complex- 
conjugated Euclidean gamma matrices. 
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[The precise meaning of this notation can be seen from the example in (|A19|) .] They are therefore only invariant under 
certain combined spin and taste rotations, and break the full spin-taste symmetry down to the discrete subgroup 
respected by the lattice theory, X SW4, 30]. In other words, these "FF(B)" operators break the continuum 
symmetry maximally. 

One advantage of this basis is that the U(1)a symmetry of the lattice theory (which is exact in the chiral limit) is 
manifest. Under this symmetry 

Q — >exp(ia[7 5 ®£ 5 ])Q, Q — ► Q exp(ia[7 5 <g> £5]) , (5) 

and each of the operators in © and Q is separately invariant. Another advantage is that there are fewer operators 
to consider; Ref. [12] had an additional six FF(A) and two FF(B) operators. 

Thus far we have only considered the single flavor theory. We now recall the generalization to n > 1 flavors [l3| . 
The quark field picks up a flavor index, Q — > Qi, and bilinear operators generalize straightforwardly, e.g. mQQ — > 
'Yl,i' m iQiQi- The generalization of the four-fermion operators is more subtle. Since these operators are the result 
of hard gluon exchange, and gluons do not change flavor, they can be written in a form containing flavor diagonal 
bilinears, i.e. 

J2Qi(is®ZT)QiQ j (is®£T)Q j . (6) 

hi 

Of course we could Fierz transform such an operator to bring it into a "flavor-mixed" form, but it is advantageous 
not to do so. This is because, in the unmixed basis 10, each bilinear must have odd distance in order to be invariant 
under the lattice axial symmetry, as discussed in the following paragraph. Thus the spin-taste structure of the FF(A) 
and FF{B) four-fermion operators remains as in and |4"|. respectively, but the flavor structure is now given by 
©. As for the single flavor case, there are two possible color contractions of each operator, but we do not show these 
explicitly as they play no role when we map the operators into the chiral Lagrangian. 

We emphasize that the result that only odd-distance bilinears appear in the flavor unmixed basis follows from 
symmetries alone. The discussion of gluon exchanges is given only to provide motivation for considering the flavor 
unmixed basis, but is not strictly necessary. Since the four-fermion operators remain in the chiral limit, they must be 
invariant under the symmetries of the lattice theory in that limit. These include U(n) axial transformations, which 
act as in except that a is now an Hermitian nx n flavor matrix, and Q a flavor vector. 4 In the lattice theory, this 
corresponds to rotating even and odd sites by opposite U(n) transformations. It is easy to check that the operators 
in J3J and Q are invariant under U (n) transformations if they are composed of flavor diagonal bilinears. On the 
other hand, compound operators with spin-taste structure such as [S x S] — [P x P] are not invariant for more than 
one flavor. We stress that it is not sufficient to consider only U(l) n transformations (in which each flavor is rotated 
independently), since this does not rule out compound operators of the form 

Yl [5i(7i ® Zi)QiQi(n ® £r)Q< - Q t hs ® ^)QiQi(i5 ® ■ (7) 

i 

These are only ruled out by their lack of invariance under general U(n) transformations. 

Finally, we comment briefly on operators of yet higher dimension. Since we aim to determine operators in S%PT 
proportional to a 4 , one might have thought that we need to enumerate all dimension 8 operators in the continuum 
effective theory (dimension 7 operators are forbidden by the axial symmetry). Fortunately this is not the case. 
The naive argument for this is that the dimension 6 operators break the continuum symmetries down to the lattice 
subgroup, and so dimension 8 operators break no further symmetries. Since it is the symmetries that are used to map 
to the chiral Lagrangian, the dimension 8 operators are not needed. This argument is, however, incomplete. What 
matters in addition is the number of taste matrices with correlated indices in the quark-level operators. If there are 
more such taste matrices, then one can build different operators in the chiral Lagrangian. We explain this point in 
Appendix lAl where in Section I A 31 we show that certain operators in the chiral Lagrangian do not arise until C(a 8 ), 
since they require four tensor taste matrices with fully correlated indices, and this does not occur unless one has an 
eight-quark operator. Now, the dimension 6 operators have at most two taste matrices, as exemplified by [A x T] in 
@. Moving to dimension 8 adds further derivatives or gluon fields, but no further bilinears, and thus no more taste 
matrices. For this reason we do not need to enumerate the dimension 8 operators. 



4 These transformations do not form a group for n > 1, but this is not important for our considerations. 
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III. CHIRAL LAGRANGIAN FOR STAGGERED QUARKS AT LEADING ORDER 



In this section we review the mapping of the continuum quark effective Lagrangian into the effective chiral La- 
grangian, working at leading order in our joint our expansion about the chiral and continuum limits. This was done 
for a single staggered flavor in Ref. , and generalized to multiple flavors in Refs. 0, 0OJ . 

If there are n flavors of staggered fermions, then, in the combined chiral-continuum limit, the theory possesses an 
SU(An)i J x SU(An)ji chiral symmetry. We assume that, as in QCD, this breaks spontaneously down to SU(An)y. 5 
This symmetry breaking pattern leads to 16n 2 — 1 pseudo-Goldstone bosons. We can collect these in the usual way 
into an SU(An) matrix: 



where $ is a traceless 4n x 4n matrix: 



S = exp(i$//) . 



( U 

IT 



TT+ K+ 

D_ K° 

K° S 



(8) 



(9) 



with 4x4 submatrices: 
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U = J2 U aTa ■ 



(10) 



and so forth. The pion decay constant, /, is normalized such that « 132 MeV. We choose to express the SU(4) 
generators in the following Hcrmitian basis: 



To, {^5: i^liVj £t} : 



(11) 



where is just the 4x4 identity matrix. It is important to retain the taste singlet meson, Ui oc tr(U), because 
with n flavors (and thus 4n tastes), only the overall SU(An) singlet decouples. Under an SU(An)L x SU(4n)j{ chiral 
symmetry transformation, E transforms as 



where L G SU{An) L and R e SU(An) R . 

We use the following power-counting scheme when determining the staggered chiral Lagrangian: 



(12) 



P 2 / A qcd ~ m/A QC D ~ a 2 A 2 QCD , (13) 
which is consistent with parameters of current simulations. The lowest order Lagrangian, which is of 0(p 2 , m, a 2 ), is 

C x = ^Tr^E^Et) - L/ 2 Tr(ME + Mtf) + a 2 V , (14) 
where M. is the quark mass matrix: 



M 



( m u I •••\ 

m d I ••• 

m s I ■■■ 

V i i 



(15) 



We do not need to worry about a possible change in the dynamics for large numbers of flavors, e.g. the loss of asymptotic freedom, 
because practical applications always involve three or less light dynamical flavors, for the which the chiral symmetry breaking pattern 
of QCD should apply. The reduction from An quarks to 2-3 sea quarks is accomplished by partial quenching, as explained at the end of 
this section. The theoretical status of symmetry breaking in partially quenched theories has been discussed in Ref. |32jl. 
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and V is the taste-breaking potential resulting from the four-fermion operators in the quark effective action. Note 
that Tr is a full An x 4n trace in both flavor and taste space. The dimensionful constant, fj,, is of 0(Aqcd), and is 
defined such that the PGB mass is given at LO by 



= 2/x- 



a 2 A F . 



(16) 



The labels i and j indicate the flavors of quarks in the PGB, here assumed to be different, while the splitting, A^, 
depends on the taste. 

The mapping of the four-fermion operators enumerated in the previous section into the mesonic operators in V is 
done by treating the taste matrices as spurions. Having worked out the implications of the symmetry under axial U(n) 
transformations at the quark level, no further subtleties arise in the mapping. The method of Ref. generalized 
to multiple flavors 0], can be used. A key result is that only the FF(A) operators, i.e. those in (JHJ), contribute to 
V, so the potential has a larger symmetry than the underlying lattice theory. We give a brief recapitulation of the 
determination of the form of V at the end of Secs. lA 2"al and lA" 2 bl as a byproduct of our extension of the methodology 
to NLO. In the notation of Ref. 113, the result is 



where 



V = U+U' 

U^Y, C ^ = CiTr^E^S 1 ) 
fc 

V 

+C 4 i^[Tr(^ ) E^ ) S) + / l .c.] 
+C 6 ^Tr^S^St) 



(17) 



(18) 



and 



C 2 y|^[Tr(^")S)Tr(^S) + h.c] 
+C 2A \ ^[Tr(^ ) S)Tr(^ ) S) + h.c] 

V 

+ c, 5^^E[ Tr ^ n)s ) Tr (^" )st )] 

V 



(19) 



(n) 

Here £j, is the 4n x 4n generalization of the taste-matrix: 








• 









• 












V ; 






■J 



(20) 



with £p the ordinary 4x4 taste matrix. This flavor diagonal form mirrors that of the underlying quark operators, 
©■ 

The potential is both rotationally and 5*0(4) taste symmetric. Therefore the PGB sector respects a larger symmetry 
than the lattice theory at 0(a 2 ). Although the 16 PGB masses (for a given choice of constituent quark flavors) are 
no longer degenerate as in the continuum, they split into five degenerate groups according to S'0(4) representations 
having tastes £5, £4, £r, £y, and £j. The taste £5 PGB is the lattice Goldstone boson and receives no mass correction 
from V because of the exact axial symmetry. Thus in (|16|l in fact only depends on the 5*0(4) representation of 
the particular taste. 
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The potential can be simplified for a single flavor using Fierz transformations. In particular, all double-trace 
operators can be transformed into operators with only single traces. This single trace basis is that used in Ref. 12]. 
Such transformations, however, are not useful for multiple flavors because they do not maintain the flavor diagonal 
form of Indeed, as realized in Ref. the double trace operators of W give rise to a previously unnoticed effect, 
namely the presence of hairpin diagrams (quark disconnected contractions) for flavor singlet mesons with non-trivial 
taste. In the LO potential, these occur only for vector and axial tastes, but, as we show below, at NLO they arise for 
all tastes. 

Given the importance of the axial U(n) transformations in restricting the form of the underlying quark operators, 
as explained in the previous section, it is worthwhile discussing these transformations at the mesonic level. They 
are simply SU{Ati)a transformations, acting on X as in O, with L = R) — exp(ia ® £5). Here, as in the previous 
section, a is an Hcrmitian n x n matrix acting in flavor space, while £5 acts in taste space. Using the commutation 
relations among the taste matrices it is easy to show that 



?5 



An) 

= 



(21) 



It then follows immediately that all the operators in the potential V are invariant under these transformations. 

We conclude by generalizing the results to partially quenched (PQ) theories, since partial quenching is often used 
in practical simulations. Following Ref. 33], one introduces valence quark flavors with ghost partners, each coming 
here in four tastes. There are then N = 4n va i + 4n sea quarks and M — 4n va j ghosts in total, where n va i and n sea are 
the number of valence and sea flavors, respectively. Note that the unquenched sea-quark sector lies within the PQ 
theory, so one loses no g enerality by considering the larger theory |3^ |. 

As noted in Ref. [Lj, the form of the chiral Lagrangian in the PQ theory is unchanged, except that the chiral 
symmetry group becomes SU(N\M)l X SU(N\M)r, traces are replaced by supertraces, and the PGB matrix is 
enlarged. For example, in a theory with two valence quarks (x and y) and three sea quarks (u, d and s), the PGB 
and mass matrices become (schematically) 



( X P+ 

P- Y 



V 



u 




K+ 


ir~ 


D 


K° 


K~ 


K° 


S 



\ 



X P+ 
P- Y J 



( m x I 



M = 



m v I 



V 



m u I 

m d I 

: m s I 
: : m x I 
: : : m y I ) 



\ 



(22) 



where each block is a 4 x 4 taste matrix. The taste matrices are also enlarged to size (N + M) 2 , but they maintain 
their flavor diagonal form, (|20|l . In fact, for the sake of brevity we drop the superscript (n) in the remaining sections. 
Finally, for practical calculations in S^PT, it may be simplest to enlarge S to be a U(N\M) matrix, and then remove 
the singlet component by adding a mass term to C x proportional to mj)St.r($) 2 , where mo is sent to infinity at the 
end |H . 



IV. CHIRAL LAGRANGIAN FOR STAGGERED QUARKS AT NLO 

We have extended the construction reviewed in the previous two sections to NLO. At this order, the new operators 
introduced by discretization are proportional to a 2 p 2 , a 2 m and a 4 . We explain in detail how we enumerate these 
operators in Appendix This requires combining appropriate numbers of taste spurions, mass spurions and Lie 
derivatives into independent operators that are invariant under the chiral group, forming parity invariant combinations, 
and reducing them to a minimal set using the equations of motion. We use the graded group theory method of Ref. j3|J 
to determine the number of linearly-independent operators. We include sources for left- and right-handed currents and 
scalar and pseudoscalar densities, so that we can calculate matrix elements of these operators. This leads to additional 
operators of 0(a 2 p 2 ) that are proportional to the sources, and thus contribute only to PGB decay constants. We 
refer to these as "source" operators. We collect all of the NLO taste- violating operators in the tables of Appendix iBl 
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In this section we focus on the pattern of symmetry-breaking exhibited by these NLO operators. Unlike the LO 

FF(A') FF(B} 

staggered potential, they come from four-fermion operators in both S 6 and S 6 . Thus they break the 50(4)- 
taste symmetry down to the lattice symmetry group. In order to understand the symmetries and symmetry-breaking, 
we examine the consequences of these operators for physical observables. We discuss the symmetries respected by 
these quantities order-by-order in the combined chiral-continuum expansion. It is useful to keep in mind the general 
form of this expansion, illustrated by the NLO PGB mass: 

K)nLO = K) LO + Ml-W + (POm* + ( Srn Z)a>m + ( S V ' ( 23 ) 

The subscripts indicate the order, in S^PT, of the terms which contribute to the physical quantity, with m indicating 
both m and p 2 . We use a similar notation for the axial and pseudoscalar decay constants, f A and / p . 6 
The LO pion mass is given in JTBJ, while the LO decay constants are 

(/ A )lo = /> (/%>=*»/• ( 24 ) 

The results of the previous section imply that the LO contributions to these quantities respect the 50(4) taste 
symmetry, since the potential V does. This has already been discussed for the PGB mass, but also holds for the decay 
constants. In fact, the LO decay constants are SU(A) symmetric, since V does not contribute to either one. The 
staggered potential also generates interesting predictions for the difference between the properties of flavor singlet 
and non-singlet particles. This difference is predicted to vanish for decay constants of all tastes, and for the masses 
of taste £5 and pions (since there are only vector and axial hairpins). 7 One of our aims is to determine at what 
order these predictions fail. 

The form of the 1-loop contributions is 

(K)i_ lKS ~ [K) LO + a2 ] 2ln K) L o (25) 

(^ P ) 1 - ioop - [K) L0 + a 2 ]lnK) L0 . (26) 

These involve only LO vertices, and thus do not break the SO (4) taste symmetry. They do, however, break the SU (4) 
symmetry of decay constants down to 50(4). Note that, unlike in continuum xPT, t ne coefficients of the logarithms 
are not proportional to the LO mass — there are additional a 2 terms from the hairpin vertices and four-pion vertices 
produced by V. The exception is in the lattice Goldstone mass, which is protected by the U(1)a symmetry. The 
one-loop contributions to the lattice Goldstone mass and axial decay constant have been calculated in Ref . [Ij, E3 ■ 
Now we turn to the analytic NLO contributions. The generic form of these contributions are as follows: 

(5f A ' P ) m >~rn, (6f A > P ) a > m ~« 2 - (28) 

(A 4 *) m2 ~m 2 +mp 2 +p\ (A^) a2m ~a 2 m + a 2 p 2 , (A 4 -) q4 ~ a 4 . (29) 

Here we begin to include the 4-pion scattering amplitude, A 4 ™ , as it is the simplest quantity showing the most general 
pattern of contributions from the NLO operators. In these expressions we use a schematic notation, leaving out 
factors of Aqcd, but distinguishing between mass and momentum-dependent terms. Note that there is no 0(a 4 ) 
contribution to f A ' P because such terms in the chiral Lagrangian contain no sources. In order to determine whether 
there are relationships between splittings, one needs to know which types of operators from Appendix IBI contribute 
to which of the corrections in I|27I29|I . This information is collected in Table [I] for single supertrace operators, and in 
Table HT1 for double supertrace operators. The former are the only contributions to the masses and decay constants of 



We use "axial decay constant" to refer to the the standard decay constant, which comes from the vacuum-to-pion matrix element of 
the axial current, (0|v4 M |7r) <x f A . In a slight abuse of terminology, "pseudoscalar decay constant" refers to the decay constant from the 
corresponding matrix element of the pseudoscalar density, (0|P|7r) <x f p . 

We are phrasing this discussion as if the hairpin vertices iterate to give only a mass shift for flavor singlets. As is well known, however, 
this is not generally true in PQ theories because valence quark loops are absent, and there are factors of 1/4 added by hand due to 
the v'Det trick. There will generally be additional, unphysical double poles in correlators 34] . Nevertheless, one can always determine 
the size of the hairpin vertices by measuring the residue of the double pole and comparing to the prediction of PQ SxPT, as discussed 
below. 
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TABLE I: Contributions of single supertrace NLO S^PT operators to physical quantities. The column labeled "5*0(4)?" 
indicates whether the contributions are consistent with S*0(4) taste symmetry. The last five columns indicate which quantities 
receive contributions from the particular operators. For further explanation see text. 
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TABLE II: Contributions of NLO SxPT operators with two supertraces to PGB masses and decay constants. Notation as in 
Table £n(5) indicates both vector and axial taste. 

flavor non-singlet mesons, while the latter gives the additional hairpin contributions for flavor singlets. The distinction 
between single and double supertraces is less significant for the scattering amplitude, so we only include it in the first 
table. 

As indicated in the tables, there is an important distinction between the contributions from underlying FF(A) and 
FF(B) four fermion operators: the former cannot break the taste S'0(4) symmetry, while the latter can. These two 
types of operator differ primarily in their index structure. Those resulting from FF(A) operators must have indices 
contracted in pairs, while those from FF(B) operators must contain more than two identical indices. This can be 
seen by comparing the 0(a 2 p 2 ) FF(A) operators in Tables IXVIII and lXVIIII to the corresponding FF(B) operators in 
Table IXIXl or the 0(a 4 ) FF(A) operators in Tables IxlClXXlTl to the corresponding FF(B) operators in Table IXXIIll 
Note, however, that the there are no 0(a 2 m) contributions from FF(B) operators because there are not enough 
indices to contract more than two at a time. Thus 0(a 2 m) operators do not break taste 5*0(4). 

An interesting feature of the pattern of 5*0(4) breaking is that the 0(a 4 ) mesonic operators resulting from double 
insertions of underlying FF(B) operators (which are listed in Table IXXlll() do not break the 50(4) symmetry of pion 
masses, despite the fact that they break 50(4) in general. This is because, when calculating tree-level masses from 
operators such as 

E Str(^£^£^£t^£t) , (30) 

two of the four E's in each operator must be replaced by the identity, so the taste matrices collapse into a form which 
is 50(4) symmetric. In contrast, one contribution to the four-pion scattering amplitude replaces each E with a pion 
field, giving an 50(4) breaking contribution to A 47T . Other contributions do not, however, break 50(4). This is 
why, in Tabled operators of type "a 4 FF(B) V are listed as both conserving and breaking 50(4), and why the latter 
contributes only to A 47r . 

This distinction between two- and four-pion contributions does not arise for the 50(4) breaking 0(a 2 p 2 ) operators. 
These violate 50(4) by correlating indices on derivatives and taste matrices, e.g. 



(31) 
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where is the covariant derivative. They produce direction-dependent contributions to masses and decay constants, 
as well as to scattering amplitudes. 

As is shown in the tables, the 0(a 2 m) operators contribute both to pion masses and to f p (since M and M> 
are sources for P), but not to f A (since they contain no derivatives). In particular, each a 2 m contribution to PGB 
properties comes with two independent mass dependencies, m x + m y (where x and y are the valence flavors) and 

m se a, The only exception is the mass of the taste £5 pion, which is guaranteed to vanish when m x + m y — > by the 
axial symmetry, and indeed has no a 2 ^ m sca contributions. Similarly, it is easy to see that it has no contributions 
proportional to a 4 . Another feature indicated in the tables is that the "a 2 source" operators contribute only to f A . 

Finally, we compare the entries in the single and double supertrace operator tables, keeping in mind that the double 
supertrace operators are responsible for splittings between flavor singlet and non-singlet PGB properties. The entries 
are almost identical. The only difference is that some of the "Y" entries in the latter are replaced by "£^(5) only" . This 
is because the a 2 m and a 2 source two supertrace terms only contain vector and axial tastes. (Recall that the same 
is true of the LO two supertrace operators in U' .) In contrast, operators of 0(a 2 p 2 ) and 0(a 4 ) with two supertraces 
contain all possible taste matrices. 

We now summarize the consequences of the previous observations, with the goal in mind to generate predictions 
of SxPT which can be tested on the lattice. We focus on the masses and decay constants, as these quantities are 
straightforward to calculate in simulations. 

1. The SU(4) symmetry in the LO decay constants of flavor non-singlet PGBs is broken down to SO(A) by the 
NLO terms resulting from FF(A) operators. 8 There is no relation between the SU(A) — > SO(A) breaking in f A 
and f p , since the a 2 p 2 and a 2 m operators contribute differently to these two decay constants. Furthermore, 
1-loop contributions to m 2 , f A and f p will give independent contributions to the splittings. Thus there are no 
predictions between splittings at this stage. 

2. The breaking of SO (A) down to T4 XI SW4 occurs first in the NLO analytic terms, and not in the 1-loop terms. 
It arises only from two types of operator — a 2 p 2 and "a 2 source" . In particular, since both and f p only 
receive S'0(4)-breaking contributions from 0(a 2 p 2 ) operators, there are relations between the splittings within 
SO (A) multiplets. The "a 2 source" operators break all such relations among splittings in f A . There are also 
predictions for the rotational symmetry breaking in the dispersion relations. These results hold separately for 
flavor singlet and non-singlet pions. 

3. The absence of taste £5 and £ M „ hairpins at LO does not hold at NLO, where there are hairpin vertices for all 
tastes. We note that the presence of a hairpin vertex for taste £5 is consistent with its Goldstone nature since 
the vertex is proportional to p 2 , and thus m 2 at the pole. 

4. The hairpin contributions to f A and f p are related to each other for tastes P and T, since there are no two 
supertrace operators with sources for either the pseudoscalar density or the axial current. This is true for both 
the S'0(4)-conserving and the SO (A)- violating contributions. 

We discuss the various relations and predictions in detail in the following section. Here we stress that none of these 
relations follows from a symmetry of the lattice theory — indeed, we have checked that all of them are broken by higher 
dimension operators in S%PT. 



In this section we work out the detailed form of the relations which follow from the particular form of the NLO 
analytic terms. 

We first study S"0(4)-taste and rotational symmetry breaking in the pion dispersion relations. This arises from the 
a 2 p 2 FF(B) operators, which are listed in Table IXIXI Of the 18 such operators, only 8 contribute to SO(A) breaking 
in single pion properties. The others, such as 



V. 



NLO RELATIONS FOR PGB MASSES AND DECAY CONSTANTS 




(32) 



We have checked that there are sufficient independent operators to completely break the symmetry down to SU(A) in this and all other 
similar cases that we mention. 
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(with p.c. indicating parity conjugate) give either vanishing, or taste symmetric, contributions to two pion properties, 
because the two pions must be drawn from the fields with derivatives acting on them. The contribution of the 5*0(4) 
breaking operators to the chiral Lagrangian is 

+Ci [Str(E9 M E t ^Ea M S t ^)+p.c] + C 13 [Str(ES^E t ^)Str(E^E t ^) + p.c] } 
+a 2 {CWStr(Ed M Ete M Std M E£ M ) + C 36A Stv(i:d^^d^ 5lM ) 

+C 41 ^Str(^Et^)Str(a M E^) + C 41A Str(9 M Et^ 5 )Str(9 /i E^)} . (33) 

The eight operators reduce to six for two pion contributions, three single supertrace and three double supertrace 
operators: 

J2 E { E - ^C 10 )Stv(d^ v d^ Ufi ) + (C 7 - 2C 13 )Str(^$^)Str(^$^)] 
+C 36V Strid^d^) + C ilv Str(c^)Str(9 M $y 

+C 36A Str(d^ 5 d^^) + C 41A Str(d M *£ M5 )Str(d M *& M )} • (34) 

We consider first the effect of these operators on the properties of flavor non-singlet mesons, and in particular those 
composed of a quark of flavor a and an antiquark of a different flavor b, where these flavors can be either valence or 
sea quarks. Examples are kaons and charged pions, both unquenched and partially quenched. 9 For these particles 
only the single supertrace operators contribute, and we find the following results for their dispersion relations: 



Ej = p 2 +m 2 I (l + 5i), (35) 

E 2 = p 2 +m 2 (l + S 5 ), (36) 

Et = (p 2 + p 2 ) + p 2 k (l + S k - 6 A ) + ml(l + 6 k ) : (37) 

E 2 = p 2 (1 + <5 4 - 5 k ) + ml(l + 5 4 ) , (38) 

E 2 k5 = (p 2 +p 2 )+p 2 k (l + 5 k5 -S i5 ) + m l 5 (l + 6 k5 ), (39) 

E 2 5 = p 2 (1 + 6 i5 - S k5 ) + m%(l + S i5 ) , (40) 

Ef m = { P f + P 2 m )(l + S lm - S k4 ) + p 2 + m^(l + S im ), (41) 

E 2 ki = [pi + P 2 n )(l + 5 ki - 6 lm ) + p 2 k + m 2 „(l + S ki ) . (42) 



Here Ep is the energy of the taste F PGB determined from the exponential fall-off of the two-point function along 
the Euclidean time direction. The quantities mp are the masses including all NLO contributions except that from 
the £0(4) breaking operators in l|34|) . They are thus labeled by indices of the 50(4) irreps: /, 5, p, p,5 and \iv. 
The additional contributions from the 50(4) breaking operators are denoted bp. Here the labels, like those on the 
energies, are those of the 8 irreps of the lattice timeslice group, which distinguishes between spatial and timelike 
indices [37| • This is the maximal splitting among tastes consistent with the lattice symmetry for states at zero spatial 
momentum. The expressions show that, at non-zero momentum, the states fall into even smaller representations, due 
to the c oup ling of the spin and taste transformations in the lattice symmetry group. These irreps have been discussed 
in Ref. Hj. 

The expressions for the £0(4) breaking mass shifts ((5m 2 /m 2 ) are 

8a 2 

Si = —-p-yCsev + C 36 a) , 
8a 2 

£5 = —p^{C 3 Qv + 36 a) , 



9 The following relations also hold for flavor neutral mesons, e.g the 7ro-like states with flavor structure (aa — 66), if flavors a and 6 are 
degenerate. 



12 



,2 



5k = -~p{C2 - 2Cio) + -p{Cmv — Czqa) , 



— — (C 2 - 2C W ) - — 2 
^ 2 8a 2 

j2 lu ' "^2 

24a 2 -~ 



<^4 = H 7^-(C2 — 2Cio) 75"(C36y — C36A) , 



<5fe5 = 75"(C2 — 2Cio) JTf(C36V — C36A) , 



+ j, 2 {Ci — 2Cio) + -p{C36V — C35A) , 
Q„2 



f 

8a 2 

7 



2 {C36V + C^a) , 



5k4, = +-^{C 3 6V + C 36 a)- (43) 



There are three 50(4) breaking splittings, those between tastes £4 and £i, between £45 and ^5, and between £^ and 
£fc4. The six independent coefficients are sufficient to give independent contributions to each of these splittings. Thus 
there are no relations predicted between these mass splittings. 

What is predicted, however, are relations between the violation of rotation symmetry in the dispersion relations 
and the 50(4)-taste breaking in the mass splittings. These can be written as 



El- 


-El 


= 1- 


p 2 +p 2 + 2p 2 


m \- 


- m 2 


(m 2 +m 2 )/2 ' 


E k5 ~ 


E l 5 


= 1- 




m t 5 ~ 


m 45 


K 5 +m 2 5 )/2' 


E fm ~ 


E l A 


= 1- 




m lm ~ 


m k4 





(44) 

In these relations mj? is the full NLO mass. Note that the relations are trivial when p = 0. The fact that the precise 
form of the mass used in the denominator on the r.h.s. does not matter at this order allows these expressions to be 
written in terms of quantitities that are directly measurable on the lattice. 

Note that neither the taste singlet nor the taste Goldstone PGBs receive rotational symmetry violating contributions 
to their dispersion relations at this order. Such terms are prohibited by the fact that these tastes transform trivially 
under rotations. Therefore they can only feel the effect of rotational symmetry breaking through operators with four 
or more derivatives [lflj . 

which are of 0{a 2 p A ) and thus contribute only at NNLO. 

Next we turn to the relations involving the decay constants, still considering only flavor non-singlet mesons. As 
explained in the previous section, the SO (4) breaking operators of 0{a 2 p 2 ) contribute to both the masses and the 
decay constants. For the pseudoscalar decay constants, as for the masses, the contribution is through wavefunc- 
tion renormalization. Since the same terms contribute to both masses and decay constants, there are three simple 
relationships between splittings: 



(ft 


-ft 




+ /f 


'/fe5 


— fib 


fP 
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f P 
Jim 


— /fc4 


fP 
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+ /£ 



1 / m 2 -m\ \ 

2 Vm 2 +m|/ 
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( fk5-Us \ = }_( 

\m + L p J 2I 



f fL- fki \ _ 1 ( m lm ~ m U \ (A(U 
^fL + ftJ 2lm 2 m + m 2 J- W 

These are likely to be the most simple predictions to test in practice. It is essential for these expressions that the 
Z-factors are S'0(4)-invariant, and therefore identical for both tastes in the expression. They can thus be tested using 
bare lattice operators, thereby avoiding O{o?) ambiguities in matching lattice and continuum operators that could 



13 



destroy the relationship. This is not true in general (e.g. for f A ) but does hold for the pseudoscalar operators 75 <g> £f- 
For example, the operators with spin-taste structure 75 (g> £j are related to 75 <S> £4 by a lattice rotation. 

This argument shows why there can be no similar relationships involving the axial decay constants, since the Z- 
factors in that case are not SO (4) invariant. For example, 745 ® £3 is a "2 link" operator while 745 (g> £4 is a "4 link" 
operator, so they are clearly not related by any lattice symmetries, and thus have different Z's. In fact, it is not even 
clear how to unambiguously calculate the quantities on the l.h.s. of eq. I|46(l for f p — ► f A , as the Z-factors are not 
universal at 0(a 2 ). Thus it is reassuring that S%PT predicts no relations involving f A , because of the "a 2 source" 
mesonic operators, which contribute SO(4) breaking contributions to f A alone. The only prediction for axial decay 
constants is a qualitative one: the splittings between 50(4) irreps should not be different in magnitude from those 
within the irreps. In other words, there should be no pattern of approximate degeneracies, unlike for the masses which 
exhibit approximate 5*0(4) symmetry. This expectation is consistent with the results of Ref. |l6l |. in which f A was 
calculated for all tastes in the quenched approximation. 

We now turn to S'0(4)-breaking predictions for PGBs for which the double supertrace, or "hairpin", operators also 
contribute. These are particles which have a flavor singlet component, i.e. uu + dd and ss in the unquenched sector, 
and xx or yy (and their ghost companions) in the valence sector. These states are mixed by hairpin contributions, 
and because, in general, they start off before mixing with different masses (due to the different masses of the quarks), 
there are no simple predictions for the properties of the resulting mixed states. The situation is yet more complicated 
if the -v^Det trick is being used. In fact, for the vector and axial tastes these complications arise at LO, due to the 
presence of hairpins, and have been addressed in detail in Ref. jl3j . 

The only theories for which the flavor non-singlet predictions still hold as written are unquenched theories with 
N > 1 degenerate flavors. In these theories there is a single flavor-singlet state, and thus no mixing. This state can 
have any taste except £/, which is, as noted above, not a PGB. The predictions for the other fifteen tastes take exactly 
the same form as those given above for the flavor non-singlets, except that there are additional contributions to the 
(5's from the double supertrace operators: 

csinglet _ 32q 2 



Ringlet = _ Ar irr_ (C7 _ 2Cl3) _ (47) 



p 

r 32a 2 
_16a 2 

T 2 



Thus the 50(4) breaking splittings for the flavor singlets are not related to those of the non-singlets, but the form of 
the predictions for dispersion relations, l|44|l . and of the relations between mass and decay constant splittings, H46[) . 
remain unchanged. 

The situation is similar if one has N degenerate flavors and uses the vT)ct trick. The only complication is that the 
disconnected quark contractions must be multiplied, by hand, by a factor of 1/4 before being added to the connected 
quark contractions. Then the predictions of 144|l and l|4tj[l apply, with the S's having the additional contributions 
given m (|4Tj> . except that N -> N/4. 

For theories with non-degenerate sea quarks, such as QCD, and for PQ theories, there are no simple relations 
involving masses and decay constants. Instead, we can obtain simple relations involving the hairpin vertices by 
considering disconnected valence correlators. These have the flavor structure 11 <-> 22, where 1 and 2 refer to valence 
flavors (x, y, . . . ). This picks out the desired hairpin contractions at the quark level. To cancel the LO decay constants, 
and to obtain a quantity with a simple dependence on Euclidean time, we take the ratio of the disconnected correlator 
to the connected (12 «-> 21) correlator. This ratio has been studied extensively in the quenched approximation. In 
the present context, it has the attractive property that the derived relationships hold even if one uses the v^Det trick. 
For simplicity, we restrict ourselves to degenerate valence quarks, i.e. m x = m 2 . 

The earlier flavor nonsinglet predictions related splittings within S*0(4) irreps in masses and pseudoscalar decay con- 
stants. As a natural generalization to flavor singlets, we therefore consider S l O(4)-breaking in disconnected correlators 
of the pseudoscalar density. For example, for tensor taste we calculate 



(Pn (0)P 22 (p = 0,t)) itm - (Pn(0)P 22 (p = 0,t)) Cfc4 
(Pii(0)P 22 (p = 0,t))^ 



(48) 



where P indicates a pseudoscalar operator with the flavor indicated by the subscript. The taste of the operators are 
shown by the subscript to the expectation value. Note that, in the denominator, any choice of tensor taste can be 
used, as the difference between the correlators for timcslicc irreps and is subleading. 
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■ X > 

LO NLO LO 

FIG. 1: NLO contribution to the {PP) flavor-disconnected correlator. The two black squares represent insertions of the 
pseudoscalar density, while the cross represents the hairpin vertex. For tensor taste, this vertex comes from both 0(a 2 p 2 ) and 
0(a 4 ) two supertrace operators, but only two of the 0(a 2 p 2 ) operators lead to 50(4) breaking. The analogous diagram, with 
pseudoscalar sources changed to axial currents, contributes to the (AA) flavor- disconnected correlator. 



We describe the calculation of l|48|l in detail, as the remaining predictions in this section all rely on some gen- 
eralization of this method. At NLO, the numerator is due to one diagram, shown in Figure ^ in which two LO 
insertions of the flavor singlet pseudoscalar density are contracted with a NLO hairpin vertex. While the full NLO 
tensor taste hairpin vertex comes from both 0(a 2 p 2 ) and C(a 4 ) operators, only two operators, CV and C13 in i|33fl . 
produce S'0(4)-breaking. 10 Thus, before Fourier transformation, the numerator of (|48|l is 

(Pn(0)P 22 (p)) 6m - (Pii(0)P 22 (p)) 5fc4 = lQa 2 fx 2 p 2 (Cr-2C 13 ) ( 1 ) , (49) 

\P ~i~ m fiu / 



1 

P z + m z p " \p z + m? 



16a V (CV - 20 3 ) { -T-—T - m% -^—^ ) \ ■ (50) 



The mass m^v is that of the non-singlet taste tensor meson with flavor 12. In this equation it is the LO mass, which 
is the same for all components of the SO(A) irrep. At NLO accuracy, it can, however, be replaced by the NLO (or 
all orders) mass of any member of the irrep. In the second line we have shown explicitly the single and double-pole 
contributions. In an unquenched theory, the former would give rise to a correction to the pseudoscalar decay constant, 
the latter to the PGB mass. The SvPT prediction is that they have a common coefficient since they result from the 
same operators. 

To denominator of the ratio Ij48(l is just the LO contribution to the connected (PP) correlator. Before Fourier 
transformation it is 

(Pn(0)P 22 (p))^ - M 2 / 2 * 2 ■ ( 51 ) 

P 2 + m 2 ^ 

We now set p = 0, Fourier transform to Euclidean time, and form the ratio, yielding: 
(Pn(0)P 22 (p= 0,t)> 6ra - (Pn(0)P 22 (p = 0,t))t» 



(Pn(0)P 22 (p = (M))^ 



C nslet (l-[l + rrv |*| ]/2) 



-C Slet (l-"Vl*l)/2, (52) 



where t is Euclidean time, and <5^4 Slct is given in H47[) with N = 1. On the first line, the 1 on the r.h.s. is from the 
single pole in the numerator, while the [1 + m F \t\ ]/2 is from the double pole. We imagine that the mass is 
first determined from the connected correlator. The prediction of S%PT is then that the constant term and the term 
linear in m^„|i| should have opposite coefficients. That it is possible to separate these two terms has been shown in 
quenched simulations [3j| . Note that for this prediction to hold one must be at large enough Euclidean times that the 
PGB contributions dominate the correlators in both the numerator and denominator of the ratio, but small enough 
that |£fc4 lglet ?TVi / i| <C 1. If the latter condition does not hold, further terms in the iteration of the NLO hairpin must 
be included. 

As noted above, the prediction 152|l holds irrespective of whether the v^Det trick has been used. This is because 
the hairpin correlators do not, at this order, access the sea sector of the theory. The only caveat is that one must 
ensure that the mass of the tensor taste valence mesons differs significantly from that of any of the tensor taste mesons 
composed of sea quarks. If not, mixing of valence and sea sectors, although of higher than NLO, can be enhanced by 
the proximity to intermediate poles. 

Another feature of the flavor-singlet prediction common with those earlier is that one can use the bare lattice 
pseudoscalar densities, since Z factors cancel in the ratio. One might be concerned that the Z factors in the numerator 



Recall that 0(a 4 ) operators are SO (4)-invariant for 2-PGB properties, and only two-supertrace operators produce hairpin vertices. 
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and denominator might differ, since the calculation of the former involves additional quark disconnected diagrams 
with intermediate gluons. One can show, however, that the difference between the Z factors vanishes in the continuum 
limit. This is actually true irrespective of the spin of the bilinear and follows because the operator has non-trivial 
taste. Thus the difference between Z's in the numerator and denominator is of 0(a 2 ), and does not contribute at 
NLO. 

We have not been able to find a prediction as simple as Q52JI for the vector or axial taste PGBs. The analogous 
50(4) breaking quantity to consider for, say, vector taste is: 

(Pn(0)P 22 (p = Q,t)) u - (P n (0)P22(f =Q,t)) 6k _ 



(Pii(0)P 22 (p = 0,i)) 



Using the methods of Refs. 0, one can determine the form of this ratio. The calculation is complicated by the 
presence of the LO hairpin vertex, which means that there are not only terms proportional to 5 B ^ ns e , but also terms 
proportional to the non-singlet £0(4) breaking quantity (64 — 5k). Furthermore, the fact that the LO hairpins must 
be iterated to all orders brings in the sea sector through intermediate propagators, and the momentum dependence 
is quite complicated. Thus we do not give explicit expressions. 

Finally, we discuss the fourth point of our summary in the previous section. We can relate the hairpin contributions 
to f p and f A for pseudoscalar and tensor taste because they only arise from operators of 0(a 2 p 2 ) and C(a 4 ). In 
particular, there are no two supertrace 0(a 2 m) operators, which would contribute to f p , or 0(a 2 source) operators, 
which would only modify f A ; all NLO hairpin operators contribute to either both or neither quantity. It is a 
straightforward exercise to calculate the contribution of these operators to f p and f A , and, when we do, we find a 
simple relation for both tastes: 



fp \ f A 1 ' \ f p 1 \ f a , 

J / hairpin \ J / hairpin \ <* / hairpin \ J / hairpin 



(54) 



The contributions of the NLO two supertrace operators are equal and opposite. The minus sign comes from the fact 
that the 0(a 2 p 2 ) terms contribute to f p only through wavefunction renormalization, but to f A also through the axial 
sources contained in the covariant derivative. Note that these relationships do involve both 50(4) conserving and 
breaking NLO contributions — in this sense they are more general than those considered above. 

These expressions can only be evaluated as written in unquenched theories with degenerate quarks. One needs to 
determine the decay constants of flavor singlet and non-singlet PGBs, and take the difference — this is 5 /hairpin- For 
non-degenerate sea quarks and/or PQ theories, however, this will not work. As above, one can instead isolate the 
hairpin contributions using disconnected valence correlation functions. We consider the following ratios of pseudoscalar 
densities and axial currents: 



(Pll(0)P22(p = 0,*))fr (^i(0)^ 2 (P = 0,*))fr 



(P 12 (0)P 2i (p = 0,t)) fT ' (Aj 2 (0)A* 1 (p = 0,t))i 



(55) 



The superscript on the axial currents indicates the Lorentz index, and not taste. Again we assume equal valence 
quark masses, mi = m 2 , and that the valence PGBs have masses differing from those of all PGBs composed of sea 
quarks. In both these ratios bare operators can be used as Z factors cancel as above. 

First consider taste £5. In this case 0(a 4 ) hairpin vertices are prohibited because of the lattice axial symmetry, so 
that the numerators of the ratios only involve two operators: 



a 2 



{C 7P Str(£d M £t£ s )Str(£td M ££ 5 ) + Ci 3P [Str(£d /i £t£ 5 )Str(£^£t£ 5 ) +p.c.]} (56) 



The contribution to the pseudoscalar ratio can be worked out as above. Only the diagram of Fig. ^ contributes to the 
numerator, and one finds at NLO 

(Pn(0)P 22 (p = 0,^ 5 = _ 16oV {C7p _ 2Ci3p) (l-m 6 \t\\ (g?) 



(Pl2(0)P2l(p = (M))fo P 



In fact, without even relating this to the (AA) ratio, we have another prediction of S^PT: the constant and m^\t\ 
terms should have opposite signs. 

For the axial current correlator, the numerator receives contributions not only from wave function renormalization 
(Fig.[TJ, but also through direct renormalization of the axial current (shown in Pig. EJ) - The flavor-disconnected axial 
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LO NLO 

FIG. 2: NLO contribution to the (AA) flavor-disconnected correlator from axial current renormalization. The single black 
square represents an insertion of the LO axial current, while the double box represents an insertion of the NLO axial current. 
It is shown as two boxes because the hairpin vertex is in the current itself. 



correlator for taste £5 is: 

P 2 I 



(Af^A^ip))^ = -16a 2 (C 7P - 2C 13P )pl j- 



Cp 2 +m2)2 2 (p 2 + m!)/' (58) 

where the first term is from wavcfunction renormalization, and the second is from axial current renormalization. If 
we set p = 0, this can be simplified to: 

(A 4 n (0)Ai 2 (p= 0,p4)> & = l6a 2 (C 7P - 2C 13P ) {2 - ™* 2 j . (59) 

Thus, in the axial correlator, the fact that only 0(a 2 p 2 ) operators contribute for taste £5 results in the absence of a 
single pole term. To test this in practice one Fourier transforms and takes the ratio to the connected correlator: 

(Ai 1 (0)Ai 2 (p = 0,t))t s 16oV ,f -l-m 5 \t\ \ 

(AU0)AUp = 0,t)U 5 ~ ~ — {C7P ~ 2Cl3P) { 2 ) • (60) 

The absence of a single pole translates into the predicted (1 + to 5 |£|) dependence. 

Now we can return to the relation between hairpin contributions to f p and f A . This shows up in the fact that 
the overall coefficients in the (PP) and (AA) ratios are the same. The prediction in l|54|) becomes the result that the 
constant terms have opposite signs. On the other hand, the linear terms must have the same sign as they correpond 
to a mass shift which does not depend on the external operators. One way to test the prediction of opposite constant 
terms is to add the two ratios: 



(P u (0)P 22 (i?=0,% 6 , (Ai 1 (0)A 2 < 2 (p = 0,t)) ( 



(Pn(0)P 21 (p = 0,t)) & (A\ 2 {Q)A\ 1 {p = 0,t)), 



OC 7715 \t\ . (61) 



The S^PT prediction at NLO is then that there is no constant term in this quantity. 

Because the tensor taste PGB is not a lattice Goldstone pion, it has 0(a i ) hairpin vertices in addition to those of 
0(a 2 p 2 ). Because of this, neither of the predictions l(57|) or l|6Tj|) apply — the constant and linear terms in \t\ within each 
individual ratio are not related. Since the C(a 4 ) hairpins lead only to mass renormalization, one expects, however, 
for the constant terms in the (PP) and (AA) ratios to be opposite, as for taste £5. This turns out to be the case. 
Thus we find 



(Pii(0)P 22 (p = 0,*))^ (Af 1 (0)A 4 22 (p = 0,t)) i 



(Pi 2 (0)P 21 (p= 0,t))^ (A\ 2 {Q)A\ 1 {p = 0,t))i 



cx rrifiu |i| . (62) 



The absence of a constant term in this quantity is really two independent predictions, one for taste £4^ and another 
for & m . 

Since the predictions l|57|). H6U|) and H61fl do not require SO(A) breaking, there can be and in fact are one-loop 
contributions to these ratios involving LO vertices. Note that, as discussed in Ref. [32], quark disconnected correlators 
receive contributions at one-loop from single supertrace vertices. We have evaluated this contributions, and find that 
they are consistent with all three predictions. Thus these are full NLO predictions of SxPT. 



VI. CONCLUSION 



Staggered fermion simulations are currently able to reach much lower dynamical quark masses than other fermion 
discretizations. However, taste-symmetry breaking is numerically significant at current lattice spacings. Thus com- 
bined chiral and continuum extrapolations incorporating discretization errors are crucial for correct extrapolation of 
physical quantities. Moreover, the \/Det trick, which must be in practice be used before taking the continuum limit, 
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may or may not change the universality class of the theory. Thus it is unclear whether one is actually studying QCD. 
We have addressed both of these issues in this paper. 

By enumerating all of the NLO operators in the staggered chiral Lagrangian, including source terms, we allow a full 
NLO calculation, including analytic terms, of PGB properties. Because 50(4)-taste symmetry breaking first enters 
at this order, these operators correctly reflect the true symmetry group of the underlying lattice action. They are 
therefore necessary, combined with one-loop contributions, for accurate extrapolations and precise determinations of 
physical quantities. 

As we noted in the Introduction, the large number of operators that appear at NLO implies that one must calculate 
a correspondingly large number of physical quantities in order to obtain predictions. Thus, for most new quantities 
one calculates (e.g. the decay constants of PGBs of different tastes), the NLO analytic terms simply give to an 
independent contribution proportional to a 2 . The most important NLO terms are then the one-loop contributions, 
which lead to non-analytic dependence on the quark mass and lattice spacing. 

The exceptions to this observation are quantities which are only non- vanishing because of SO ( A) breaking. These 
only recieve contributions from the NLO operators we have enumerated, and, in fact, from only a handful of these. 
There are no loop contributions at NLO. This allows us to make a number of testable predictions relating SO (4)- 
breaking splitting in PGB masses, decay constants, and dispersion relations. These predictions can be tested by 
calculating only two-point correlators of unrenormalized lattice operators. For the hairpin operators we suggest a 
method similar to that used successfully in quenched simulations. Thus we hope that it will be practical to carry out 
these tests. 

We also find a few predictions for quantities that do not involve SO(4) breaking, but reflect the structure of S%PT. 
These are NLO analogues of the absence of the taste tensor hairpin at LO. We stress that these predictions, like 
those above, are not consequences of lattice symmetries, and indeed are violated at NNLO or higher order in the 
chiral-continuum expansion. In a similar vein, we note that not all operators in the effective chiral theory which are 
consistent with the lattice symmetries and power counting actually appear at a given order. One must enumerate 
quark-level operators and then match these onto mesonic operators, and this leads, at NLO, to restrictions on the 
contraction of indices. 

In our view, the most important applic ation of the predictions we have given is to test the applicability of the 
effective chiral theory when one uses the v^Det trick, and by so doing to test the assumption that this t rick does not 
modify the continuum limit. Most of our predictions arc unchanged in form when one implements the vTJet trick by 
hand in the effective theory. This is both good and bad. It is good because it makes the predictions somewhat more 
theoretically robust; it is bad because one would like to calculate quantities which directly check the factors of 1/4 
put in by hand. In any case, the key point is that if SxPT fails to describe r esult s from numerical simulations with 
small quark masses, this would shed serious doubt on the correctness of the v^Det trick. Testing our predictions will 
thus provide further empirical evidence either for or against the validity of this trick. 
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APPENDIX A: DETERMINATION OF 0(a 2 p 2 , a 2 m, a 4 ) OPERATORS 

Here we determine all operators in staggered chiral perturbation theory which are NLO in our power counting and 
arise from discretization errors. We do not discuss the remaining NLO operators, which are of 0(p 4 ,p 2 m, m 2 ), since 
they are unchanged from continuum chiral perturbation theory. 

We first discuss the chiral operators of 0(a 2 p 2 ,a 2 m). These arise from a single insertion of the quark- level 0(a 2 ) 
operators. Then we turn to the C(a 4 ) chiral operators, which require either two insertions of quark-level 0(a 2 ) 
operators, or a single insertion of an 0(a 4 ) operator. We close this subsection with a brief note on operators which 
do not arise until higher than NLO in our expansion. Finally, we consider the 0{a 2 p 2 ) operators involving covariant 
derivatives acting on taste spurions, which lead to additional contributions to the vector and axial currents. 



1. Single insertion of gluonic operators and fermion bilinears 

Most of the dimension six gluonic operators and fermion bilinears do not break continuum symmetries, and can 
only generate corrections proportional to terms in the continuum chiral Lagrangian. Thus at NLO they only give rise 
to the two operators in the LO chiral Lagrangian, 

Str(<9, I £G\ 1 I] t ) and Str(M f S + MS f ) , (Al) 

multiplied by independent coefficients of size a 2 . In other words, they generate independent 0{a 2 ) corrections to 
the parameters / and /i. As we will see, many additional contributions to the two operators in (|A1(I arise from the 
four-fermion operators discussed in the following sections. However, these, too, can simply be absorbed into the 
unknown coefficients multiplying each operator. 

Here we use mass spurions which transform in the same manner as the £ field under a chiral symmetry transfor- 
mation: 

M -> LMPJ , Aft -> RM^L ] . (A2) 

These spurions also function as sources for the scalar and pseudoscalar densities, with M = s + ip, and s and p general 
Hermitian fields. In the absence of external fields we set s — ► M. and p — ► 0. 

The remaining two gluonic and bilinear operators violate Euclidean rotation symmetry, 

J2^(D^DM, ^Q( 7m ®1)££Q, (A3) 

although they are taste symmetric. They can give rise to rotationally non-invariant mesonic operators, but this 
requires four derivatives |l2| . e.g. 

Str(d^d^d^d^) . (A4) 
These operators are of 0(a 2 p 4 ), and thus NNLO, which is one higher order than we consider here. 



2. Single insertions of operators from S^ F ' A ' 

FF( A) 

It is important to remember that the four fermion operators in S 6 are invariant under rotations and 5*0(4) 
taste transformations, so they must map onto mesonic operators that also respect these symmetries. For discussion 
purposes, we divide them according to how they transform under SU(4N\4M)l x SU{AN\AM)r chiral rotations. 



a. Operators with spin structure V or A 



Here we consider mesonic operators arising from single insertions of four-fermion operators of the form [V x F] and 
[Ax F], where F = S, P,T indicates the taste of the bilinears. The chiral structure of these operators is 

O f = ±J2 (Q^(7m ® Fr)Q r ± Ql{ la (8 F L )Q L ) 2 , (A5) 
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1 x \7\7\ + 1 



FIG. 3: The product of two bifundamentals in SU(4N\4:M). The result also applies for SU(N) for N > 3 (with the dashed 
lines removed). 



1 x 



Symm. 



\7^m + 



FIG. 4: The symmetric product of two bifundamentals in SU(4N\4M). The result also applies to SU(N) for N > 3 (with the 
dashed lines removed). 



where and the upper and lower signs correspond to spins V and A, respectively, and Fr.l are Hermitian taste 
matrices. To determine the corresponding mesonic operators, we promote the taste matrices to spurion fields, with 
transformations chosen so that Of is invariant under chiral transformations: 

F L ^LF L L\ F r ^RF r E). (A6) 

At the end we will set Fl = Fr = F, where F is the specific taste matrix appearing in the operator. 

We first construct the resulting 0(a 2 p 2 ) operators. These must be chiral singlets, respect parity (Fl <-> Fr and 
£ £ T ), be quadratic in F (odd powers of F are forbidden by the F — ► — F symmetry and quartic terms are higher 

FF( A) 

order in a), and contain two derivatives, S 6 is rotationally invariant, so the indices of the derivatives must be 
contracted with each other, and cannot be correlated with any of the indices associated with the taste matrices. It is 
easiest to construct linearly-independent, chiral singlet operators using elements with the same chiral transformation 
properties. We use left-handed building blocks, the choices for which are: 

F L -> L(F L )tf 

E0„Et ^ £(Ea M Et)Lt. (A7) 

We can now use the graded group theory method of Ref. 36] to determine the number of linearly independent 
0(a 2 p 2 ) operators in a general SU(4N\4M) partially quenched theory. Because all of the operators in (|A7I) transform 
as bifundamentals under the left-handed chiral group they all potentially have both adjoint and singlet components. 
An operator with two taste spurions therefore comes from the product of an adjoint plus singlet times an adjoint plus 
singlet, as shown in Fig. fusing the graded Young tableaux notation of Refs. If the two taste spurions are 

identical (both Fl or Ei^E 7 ), it will only come from the symmetric part of this product, shown in Fig. EI The left- 
handed Lie derivative has vanishing supertrace, however, and therefore no singlet component. An operator with two 
Lie derivatives, which necessarily have their indices contracted and are therefore identical, comes from the symmetric 
product of two adjoints, shown in Fig. El Thus the number of linearly-independent operators is equivalent to the 
number of singlet representations contained in the product of Fig- Ej which comes from the Lie derivatives, and either 
Fig- El °r El which comes from the taste spurions. 

Consider first operators composed of two Lie Derivatives, one Fl, and one EFrEL This requires counting singlets 
in the product of the representations in Figs. El and El Since all of the representations in Fig. El are self-conjugate, 
each can form a singlet with the identical representation in Fig. El This leads to eight linearly-independent operators. 
We choose the operator basis shown in Table ITTTI in which operators have been simplified using the anti-Hermiticity 
of the Lie derivative and EE' = 1. Each listed operator has the same coefficient with the same sign for both V and 
A spins because the overall sign and the cross-term sign in i|A5|) cancel. Since group theory only considers chiral 



Symm. 

FIG. 5: The symmetric product of two adjoints in SU(4N\4M). The undotted boxes are fundamental representations and 
the dotted boxes are anti-fundamental representations. The result also applies to SU(N) for N > 3 (with the dashed lines 
removed) . 



f - + EH + 1 
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Operator 



Keep'i 



Str(^Et^EEtFiEFfl) 
Str(a M E t F i a M EFfl) 
Str(a M E9 M E t EF ii E t Fz,) 
Str(a M E9 M E t F i )Str(F fl ) 
Str(a M E t a M EF JJ )Str(F L ) 
Str(a M E t a M E)Str(F i EF fl E t ) 

7. Str(E9 M E t F L )Str(E t a M EF fl ) 

8. Str(a M E t a M E)Str(F i )Str(F B ) 



Yes - combined with 3 
Yes 

Yes - combined with 1 
No 
No 
Yes 
Yes 
No 



TABLE III: The eight linearly- independent 0(a 2 p 2 ) operators with two derivatives and two different taste spurions. Note that 
two of them must be combined and some can be neglected because of other considerations, as discussed in the text. 





Operator 


Keep? 


9. 


±Str(9 M Ea M E t F L F i ) ± Sir (d ^ d„ZF R F R ) 


No 


10. 


±Str(E9 fl E t F L Ea M E t F L ) ± Str(E t a M Ef fl E t 9 M EF H ) 


Yes 


11. 


±Str(9 M E t a M E)Str( J F L F L ) ± Str^Et d M E)Str (FrFh) 


No 


12. 


±Str(9 M E9 fl E t F L )Str(F L ) ± Str(cVE%EFfl)Str(F fl ) 


No 


13. 


±Str(Ea fl E t F i )Str(Ea M E t F i ) ± Str(E t 9 M EF fl )Str(E t 9 M EF fl ) 


Yes 


14. 


±Str(^E t a M E)Str(F L )Str(F L ) ± Str(a, l E t 9 fl E)Str(F fl )Str(F J? ) 


No 



TABLE IV: The six linearly-independent 0(a 2 p 2 ) operators with two derivatives and two identical taste spurions. The upper 
signs correspond to spin V, while the lower signs correspond to spin A. Some operators can be neglected because of other 
considerations, as discussed in the text. 



transformation properties, we must impose other symmetries by hand. The operators 1 & 3 and 4 & 5 transform into 
each other under parity, so we include their sums, each with a single undetermined coefficient, in the staggered chiral 
Lagrangian. Operators 4, 5, and 8 are all proportional to Str(F), so they vanish unless F = £/. When this is the case, 
however, they are proportional to the LO kinetic term, which we have already included in the previous sub-section. 

Similarly, six linearly-independent operators can be made out of two Lie Derivatives and two Fl's or Fr's. These 
are shown in Table IIVI Here the minus signs in l|A5Jl do not cancel, and the coefficients of operators with V and A 
spins have the same magnitude, but opposite sign. Operator 12 vanishes unless F = £/, and then is proportional to 
the LO kinetic term, as are operators 9, 11, and 14 for all tastes, so we can neglect them. After combining both lists, 
we are left with six new linearly independent 0(a 2 p 2 ) operators. 

We note in passing that the fourteen operators in Tables ITTT1 and Hvl are those obtained simply by building operators 
out of the appropriate elements and using the cyclicity of the supertrace. There are no relationships between these 
operators in the PQ theory, unlike for chiral 5/7(2) and SU(3) theories. Typically, such relationships are found 
using Cayley-Hamilton relations, as in Refs. ^(j an d ^U- There is, however, no graded analog of Cayley-Hamilton 
relations, since the superdeterminant is not a finite polynomial |39| . 

To obtain the actual operators in the chiral Lagrangian we set Ft, = Fr = F, a specific taste matrix, here S, P or T. 
Since rotational and SO{A) taste symmetry are unbroken, the indices on the derivatives must be contracted together, 
as must those on the taste matrices. With these rules, it is easy to see that the specific 0(a 2 p 2 ) operators in the Sx>C 
generated by operators in 5 ,FF ( A ) with spins V and A are those shown in Table IXVlII We have dropped the ± signs 
at this stage because we do not know the relative size of the coefficients of the underlying four-fermion operators with 
spins V and A, and therefore cannot make use of the fact that their mappings into the chiral Lagrangian are related. 
We have also dropped the operators resulting from taste S, since they are either of the same form as the LO chiral 
Lagrangian, or vanish because Str(£<9 M £t) = 0. 

We determine the 0(a 2 m) operators using the same methodology. We use the mass spurions introduced above, 
which lead to two new left-handed objects: 

MS t -> LAfE'L* , SM t -> LEM+L* . (A8) 

Like the taste spurions, these objects are also bifundamentals under the left-handed chiral group. Operators of order 
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Operator 




Keep? 




15. 




Yes 


- combined 


with 22 


16. 


Str(FLEAF EFi?E T ) 


Yes 


- combined 


with 21 


17. 


ca i f t~ i v"! n t\rij / 7i /ft \ 

Str^FLEFflE 1 )Str(EM 1 ) 


Yes 


- combined 


with 23 


18. 


Str(F L EM t )Str(F fl ) 




No 




19. 


Str(F L )Str(F H M t E) 




No 




20. 


Str(FL)Str(FH)Str(EM f ) 




No 




21. 


Str(FjjE T ME T FLE) 


Yes 


- combined 


with 16 


22. 


Str(F fl E T F L M) 


Yes 


- combined 


with 15 


23. 


StrfFV y!F 3 yi'h' ; !trl'A//"y:'M 


Yes 


— pom riinpH 


with 17 


24. 


Str(F L ME t )Str(F fl ) 




No 




25. 


Str(F L )Str(F H E t M) 




No 




26. 


Str(Fz,)Str(F H )Str(ME t ) 




No 





TABLE V: The twelve linearly-independent 0(a 2 m) operators one mass spurion and two different taste spurions. Many can 
be neglected, and the rest must be combined because of symmetry considerations. 





Operator 


Keep? 


27. 


±Str(F i F i EM t ) ± Str(F B F iJ E t M) 


No 


28. 


±Str(F L F L )Str(EM t ) ± Str(F ii F J j)Str(E t M) 


No 


29. 


±Str(F z ,)Str(F L EAf t ) ± Str(F ii )Str(F fl E t M) 


No 


30. 


±Str(F L )Str(F I ,)Str(EAf t ) ± Str(F ii )Str(F ii )Str(E t M) 


No 


31. 


±Str(F z ,F L A/E t ) ± Str(F«F«Af +E) 


No 


32. 


±Str(F L F L )Str(A/E t ) ± Str(F ii F i? )Str(M t E) 


No 


33. 


±Str(F L )Str(F I ,AfE t ) ± Str(F ii )Str(FjjM t E) 


No 


34. 


±Str(Fz,)Str(F L )Str(AFE t ) ± Str(F fl )Str(F J? )Str(M t E) 


No 



TABLE VI: The eight linearly-independent 0(a?m) operators with one mass spurion and two identical taste spurions. The 
upper signs correspond to spin V , while the lower signs correspond to spin A. All of these operators can be neglected because 
of additional considerations. 

(D(a 2 m) will contain two taste spurions and one mass spurion, and be chiral singlets that respect parity. They will 
therefore come from the singlet representations contained in the product of three bifundamentals. Such operators 
may contain either two identical or two different taste spurions. 

If the taste spurions arc different, the operators simply come from the product of two bifundamentals, Figure |3 
with a third, which clearly contains six singlets. However, because there are two types of mass spurions, there are 
twelve corresponding operators in Table Ivl Each of these operators must be combined with its parity conjugate to 
form a parity-invariant operator, so there are in fact only six operators in the chiral Lagrangian: (15 + 22), (16 + 
21), (17 + 23), (18 + 25), (19 + 24), and (20 + 26). Three of these operators are proportional to Str(F), and so 
vanish except for taste S, but then reduce to the LO operator Str(£Aff + M£l). Thus we are left with only three 
new 0(a 2 m) operators. 

In fact, only two of these three are new if we consider S-matrix elements. Operator (15 + 22) reduces to the LO 
operator if M is set to the mass matrix, since M. then commutes with Ft, and Fr. The operator does, however, 
contribute SU (4) taste breaking contributions to pseudoscalar matrix elements, for then M is a source and need not 
be proportional to the identity in taste space. 

If the operators contain identical taste spurions, two of the three bifundamentals must be symmetrized, as shown 
in Figure |3J The overall product contains four singlets, and therefore generates the eight corresponding operators in 
Table IVTI However, using the fact that F 2 — 1 for all tastes as well as Str(F) = except for £/, it is easy to see that 
all of these operators reduce to Str(SMT + MET) and can be neglected. 

Turning these generic structures into operators in the staggered chiral Lagrangian is now straightforward: simply 
let Fl = Fr = F. We do not set M = M> = A4, however, so as to allow the derivation of (pseudo)scalar matrix 
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elements. The linearly-independent 0(a 2 m) operators corresponding to 5 6 operators with spins V and A are 
shown in Table IXVIII 

Finally, we briefly recall the construction of the LO operators of 0(a 2 ), so as to make contact with the discussion 
in the main text. These operators contain two taste spurions, and the only non-trivial SU(N\M) singlet is 

Str(fiEF fl Et). (A9) 

This leads to the operators 0\ and Oq in the potential U in (|18|) . 

b. Operators with spin structure S or P 

Here we consider single insertions of operators [S x F] and [P x F], where the taste can be F — V or A. The chiral 
structure of these operators is 

0' F = (0l(l® Fl)Q r± 0^(1 ®Fr)Ql) 2 , (A10) 

where and the upper and lower signs correspond to spin S and P, respectively, and Fl,r are taste matrices. To 
determine the corresponding mesonic operators in the SxC, we promote the taste matrices to taste spurion fields 
transforming as 

F L -> LF L R\ Fr -» RFrL^ . (All) 

Note that these transformations do not maintain the Hermiticity of the original taste matrices. We can, however, 
consistently impose the relation F^ — Fr, and so do not need to consider F^ R as additional variables. The left-handed 
basis for constructing the mesonic operators becomes 

Fl^ , EFr, E9 M E t , M£t, SJlft. (A12) 

The counting and construction of operators with the correct taste structure is identical to the [V, A x F] case, 
although the taste spurions are different. The primary technical difference arises when we set Fl = Fr = F, because 
the relative minus signs within compound operators are different. 

As before, there are fourteen 0(a 2 p 2 ) operators which can be constructed from two taste spurions and two Lie 
derivatives. These are shown in Table IVTT1 Notice that the ± factors differ from those for [V,A x F] because (|A10|) 
only has an internal minus sign, not an overall one. Two pairs of operators, 35 & 37 and 38 & 39, transform into 
each other under parity, so we include their sums, each with a single coefficient, in the chiral Lagrangian. Operators 
(35 + 37) and 40 both reduce to the LO kinetic operator after setting Fl = Fr = F. This leaves ten new linearly 
independent 0(a 2 p 2 ) operators. 

In addition, there are twenty 0(a 2 m) operators that contain one mass and two taste spurions. These are listed 
in Table IVIIII Six pairs of operators transform into each other under parity, and must therefore be included in the 
Lagrangian with a single coefficient: 49 & 56, 50 & 55, 51 & 57, 52 & 59, 53 & 58, and 54 & 60. However, operators 
(49 + 56), (50 + 55), and (51 + 57) reduce to the LO mass term after setting Fl = Fr = F. Thus there are eleven 
new linearly-independent, parity-invariant 0(a 2 m) operators. Three of these, (52 + 59), 61, and 63, are pure source 
terms, however, since they reduce to LO operators when M = = M.. 

To turn these generic taste structures into actual operators in the Sx£, we set Fl = Fr = F, with F being either 
V or A. However, we choose to leave M and M' as generic sources which can be set to the quark mass matrix when 
necessary. As before must contract indices on derivatives and taste matrices separately, so as not to break rotational 
or SO (4) taste symmetry. The resulting operators are listed in Tabic IXVlIII 

Finally, we recall the construction of the LO operators of 0(a 2 ). In this case there are three non-trivial SU(N\M ) 
singlets: 

Str{F L ^F L ^)+P-c, (A13) 
Str(FiEt)Str(F L Et) + p.c. , (A14) 
Str( J F L E t )Str(EF it ) . (A15) 



The first of these leads to O3 and O4 inU, (JTSJ, the second to O2V and O2A m W, (|19fl . and the third to 0$a and 
5B in W. 
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Operator Keep ? 

35. ±Str(a M Ea M E t F L F R ) ~ No 

36. ±Str(S9 M E t F i E t 9 M SF R ) Yes 

37. iStr^Et^EFRFi) No 

38. ±Str(a M E t a M EE t F L )Str(EF ii ) Yes - combined with 39 

39. ±Str(a M E9 M E t EF i{ )Str(Fi,E t ) Yes - combined with 38 

40. ±Str(a M Eta M E)Str(F i F JJ ) No 

41. ±Str(a M E t F L )Str(a M EF fl ) Yes 

42. ±Str(a M E t a M E)Str(F L E t )Str(EF ii ) Yes 

43. Str^E^EtFtEtFiE 1 ") +Str(a M E t a M EF fi EF iJ E) Yes 

44. StT(d^F L d^F L )+StT(d^T,F R d„T,F R ) Yes 

45. Str(a M E t 9 M E)Str(F L E t F L E t ) +Str(a M E t a M E)Str(EF ii EF J? .) Yes 

46. Str(a M E9 M E t F L E t )Str(F L E t ) +Str(a M E t a M EF fi E)Str(EF R ) Yes 

47. Str(a M E t F L )Str(a M E t F L ) + Str(a M EF i{ )Str(9 M EF J? .) Yes 

48. Str(a M E t 9 M E)Str(Fz / E t )Str(F i E t ) + Str(^E t a M E)Str(EF fi )Str(EF fi .) Yes 

TABLE VII: The fourteen linearly-independent 0(a 2 p 2 ) operators with two derivatives and two taste spurions. When present, 
the upper and lower signs correspond to spins S and P, respectively. 

c. Operators with spin structure T 
The spin T operators, [T x V] and [T x A], have almost the same chiral structure as those of spins S and P: 

°f = J2 (Ql(^ ® Fl)Qr) 2 + (oifr^ ® ^fl)Qi) 2 • (A16) 

The only difference is that there are no cross terms between Ft, and ffl. This means that the corresponding chiral 
operators are identical to those generated by [S 7 P x V, A], except that those coming from FlFr cross-terms are 
absent. Thus the only effect of single insertions of four-fermion operators with spin structure T is to change the 
(unknown) coefficients of some of the mesonic operators already listed. 



3. Single insertions of operators from Sg F< - B ^ 

The four-fermion operators in Sq F ^ break both rotational symmetry and the remaining SO (4) taste symmetry 

repected by operators in Sq F ^ A \ Thus they can map onto mesonic operators that also break these symmetries. 
Indeed they only map onto such operators, because they are constructed to have no taste singlet component. In this 
subsection we construct all such operators resulting from a single insertion of S^^ . 

We begin with a general comment. To break the rotational and/or SO (4) taste symmetries requires that the mesonic 
operator have more than two repeated indices. This can be accomplished at NLO either by having two derivatives and 
two taste spurions, or by having four taste spurions. The former choice leads to the 0(a 2 p 2 ) operators exemplified by 

]T Str(S^Et^St^S^) , (A17) 

and discussed in this subsection. The latter choice gives 0(a 4 ) operators such as 

J2 E Str(^£^E^St^£t) , (A i 8) 

which will be discussed separately below. It is not possible, however, to break the symmetries with two taste spurions 
and a mass spurion. Thus there are no new 0(a 2 m) operators produced by S FF< " B \ 
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Operator 

49. ±Str(F L F R SM t ) 

50. ±Str(F L A/ t SF fl ) 

51. ±Str(F L F fl )Str(EM t ) 

52. ±Str(F L Af t )Str(EF R ) 

53. ±Str(F L E t )Str(EF fi .EM 1 ') 

54. ±Str(F L E t )Str(EF fl )Str(EM t ) 

55. ±Str(F L F fl ME t ) 

56. ±Str(F L E t MF B ) 

57. ±Str(F L F fl )Str(ME t ) 

58. ±Str(F L E t ME t )Str(EF R ) 

59. ±Str(F L E t )Str(F fl M) 

60. ±Str(F L E t )Str(EF H )Str(ME t ) 

61. Str(F L E t F L A/ t ) + Str(F R ZF R M) 

62. Str(F L E t F L E t )Str(EM t ) + Str(F fi .EFflE)Str(E 1 'M) 

63. Str(F L E t )Str(F L AF t ) + Str(F fi .E)Str(F Ji M) 

64. Str(F L E t )Str(F L E t )Str(EM 1 ') + Str(F H E)Str(F R E)Str(E 1 'M) 

65. Str(F L E t F i E t ME t ) + Str(F H EF fi .EM 1 'E) 

66. Str(F L E t F L E t )Str(A/E t ) +Str(F fl EF H E)Str(M t E) 

67. Str(F L E t )Str(F L E t ME t ) + Str(F fi .E)Str(F Ji EM 1 'E) 

68. Str(F L E t )Str(F L E t )Str(ME 1 ') + Str(F fl E)Str(F R E)Str(A^ t E) 



Keep? 
No 
No 
No 

Yes - combined with 59 
Yes - combined with 58 
Yes - combined with 60 

No 

No 

No 

Yes - combined with 53 
Yes - combined with 53 
Yes - combined with 54 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 



TABLE VIII: The twenty linearly-independent 0(a 2 m) operators with one mass spurion and two taste spurions. When present, 
the upper signs correspond to spin S, while the lower signs correspond to spin P. 



Although the operators in S e appear more complicated than those in S e , because of the coupled spin and 
taste matrices, we can actually reuse the bulk of our work from the previous subsection. This is because we can still 
use spurion analysis and group theory to determine the chiral structure of the mesonic operators. We just need to 
be more careful about the index structure and the specific taste matrices. We illustrate the procedure by working 
through the mapping of [V^ x T M ] in detail; the method then can be generalized straighforwardly to the remaining 

operators in Sq F ^ . 

We begin by recalling the definition of this operator: 

Wit x T M ] = X) {^(7^ ® Z^QQiln ® &i*)Q - Q(7m ® £/«/6)Q<?(7m ® W)<9} ■ (A19) 

The second term in this expression removes the taste singlet component. It is, however, cumbersome and unnecessary 
to keep both terms, so we first simplify using 

i (\y„ x + [v x t]) = Q<M> ® MQQfri* ® MQ ■ ( A2 °) 

Since [V x T] is from S e this addition has no impact on matching to operators which break rotation and/or taste 
50(4) symmetries. The chiral structure of the operator in l|A20() is given by 

0(jjl) f = ± [Q^( 7 „ ® F(jjl) r )Q r ± Ql( 7/i (g> F{pL) L )Q L ] 2 , (A21) 

which has the same form as for [V x T], l|A5|l . except that the taste spurions depend on /i, and the sum over /j, has 
been removed for now. The mapping then proceeds as follows: first, determine the mapping of 0{h)f onto mesonic 
operators, for fixed /i; next, set the taste spurions to their appropriate values; and, finally, sum over /i and any other 
remaining indices. 

The group theory required to determine the independent mappings of <D((i)f onto mesonic operators is identical 
to that for [V x T], since the spurions F{^l)r,l must transform like Fn t L- Thus one obtains the types enumerated 
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in Tables FTTH and llVl with the provisos that -Fr,l — > F(p)r,l and that we must be careful with Lorentz (or, more 
precisely, hypercubic group) indices. In particular, we do not yet know how the indices /i in the operators in the 
tables are connected with the index /i in 0(h)f- To work this out we need to determine how 0(h)f transforms under 
rotations while keeping the spurious fixed. This may seem confusing because of the index \i in F(h)fi,Li but once, for 
a given /i, we convert the taste matrices into spurion fields, they are to be treated as rotationally-invariant (pseudo- 
scalar fields. Then, since 0(\x)f has two vector indices, it has a component which is a singlet under hypercubic 
rotations (corresponding to summing over jj,), and a component which is part of a 3-d hypercubic irrep (the diagonal 
part of the two-index symmetric tensor representation of the Euclidean rotation group) . The former is not of interest 
here, since the sum over /i decouples the spin and taste indices. It is the latter which we need to match onto mesonic 
operators. 

It is straightforward to do this matching. For each of the operators in Tables UTTl and llVl we need to replace "<9 M <9 M " 
with "dpdp — (1/4) X^p^p^p"- This projects out the singlet component. The singlet component has decoupled spin 
and taste indices, however, so we can drop it if we wish, and keep only . Doing so gives the following simple 

prescription: keep the operators in Tables ITnl and Hvl as is, but substitute Fr.l — ► F([i)r,l- In this way the spin 
and taste indices become correlated. Now that the operator mapping is done for fixed /i, we can set F (ji) to the 
appropriate taste matrix (in this case ^p V ) and then sum over the remaining indices (in this case both /i and v with 
the constraint that ^ v). For example, after completing the three steps of the mapping procedure outlined above 
for operator 1 ( Table ITTTj) we obtain 

Str^S^EE+^E^) . (A22) 

p i^#p 

This operator, like the Sg F ^ operator it came from, breaks the symmetries down to the lattice spin-taste group 

r 4 x sw 4 . 

Carrying out the same procedure for all the types in Tables ITTll and Hvl leads to the mesonic operators collected in 
the upper panel of Table IXIXI Note that the mapping of the four-fermion operator [Ap x Tp] leads to the same set 
of mesonic operators as those from [Vp x Tp) . 

The mapping of operators [Tp x Vp] and [Tp x Ap] can be done similarly. The chiral structure is given by 

0"(h)f = J2 Ql(^ ® F(h)l)Qr Qr~(^ ® F(h)r)Ql , (A23) 

where the taste spurions depend on /i, and /i is not (yet) summed. 11,12 On the other hand, v is summed, so the 
rotational properties of the operator are the same as those of 0(/j,)f discussed above. Thus when it is mapped onto 
mesonic operators with two derivatives their indices can both be set to The types of mesonic operators that arise 
are as for [S x V, A], except that here we need keep only those with one Fl and one Fr, namely operator types 35-42 
in Table IVTT1 Furthermore, here we must substitute Fr^l — ► F(h)r,l (with F(h)r.l — ¥ £p or £p5 in the end). The 
resulting specific mesonic operators are collected in the lower panel of Table IXIXI 

In Ref. [TJJ it was claimed that one of the operators generated by [T„ x V^] at 0(a 2 p 2 ) was 

Y Str [(dp - d v )^ p ] Str [(dp - 3,)^] . (A24) 
p#^#p#p 

This operator does not, however, appear in our lists. In fact, it is not invariant under the "spin-taste locked" rotation 

d v -> dp, dp -> -d u 
lv -> 7m ' 7m -» ~lv 

- £m- Cm - ■ (A25) 

which corresponds to a lattice SWi transformation Thus it does not arise. 



Note that the chiral structure of [Tp X Vp] is not the same as that of [T X V]. 
12 We thank Claude Bernard for pointing out an error in our discussion of this operator in earlier version of this paper. 
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FIG. 6: The product of two adjoints in SU(kN\kM). The result also applies to SU(N) for N > 3 (with the dashed lines 
removed) . 



4. Double insertions of four-fermion operators 

Mesonic operators of C(a 4 ) can arise from the combination of two four-fermion operators, each of which brings a 
power of a 2 , or by mapping a single 0(a 4 ) quark-level operator. However, as noted in Sec. [HJ the 0(a 4 ) operators 
do not break any additional symmetries and do not contain any more correlated taste indices, so they do not lead to 
additional mesonic operators. 

Thus we consider only the mesonic operators which result from two insertions of four-fermion operators. These two 
insertions need not have same spin or taste, and may come from either or Sg F . We first consider pairs of 

FF(A) operators, separating the various possible spin combinations. We then generalize these results to the case in 

which one or both of the operators is in S 6 . 



a. Spin V or A four-fermion operators with spin V or A four-fermion operators 

Operators of 0(a ) must contain four taste spurions - two from each four-fermion operator. The taste spurions 
corresponding to operators with spins V and A are Fl and F R , which in the end can be set to S, P or T. In fact, 
as we have seen above, the taste singlet always leads to operators which either vanish or are of the same form as LO 
operators. We have checked that this remains true for 0(a 4 ) operators, and so drop taste S from the beginning. Thus 
we can set Str(i*j,) = Str(F^) = 0, which simplifies the operator enumeration. 

There are five parity-invariant ways to combine the four taste spurions: 

{F L F R )(F' L F' R ) 
(F L F L ){F' L F' R )+p.c. 
(F L F R )(F' L F' L ) + p.c. 
(F L F L )(F L F' L )+p.c. 

(F L F L )(F R F R )+p.c. (A26) 

Here we use primes to distinguish between spurions from the separate four-fermion operators because they may have 
different tastes. To determine the number of linearly-independent chiral operators, we use graded group theory as 
before. For simplicity, we build operators out of the left-handed basis defined in (|A7|) . Both Fl and SFrSt are now 
supertraceless, so they have no singlet component and transform as adjoints under the chiral symmetry group. The 
corresponding 0(a 4 ) operators are therefore chiral singlets contained in the product of four adjoints, with various 
symmetries imposed based on the handedness of the four spurions. 

Operators of the first type, (FlFr)(F l F r ), contain four different taste spurions, and therefore possess no additional 
symmetry. The number of linearly-independent operators of this sort is simply the number of singlets in the product 
of four adjoints. (The product of two adjoints is shown in Fig. 0) There are nine such operators, and they are given 
in Table ITXI Most, however, reduce to LO operators or field-independent constants after setting F and F' to specific 
taste matrices; only three are new. 

Operators of the last type, (FlFl)(F r F r ) +p.c, contain two pairs of identical taste spurions. They therefore come 
from the singlet irreps in the product of one symmetric product of two adjoints with another symmetric product of 
two adjoints, i.e. the product of Fig. with itself. There are four such singlets; they correspond to the operators in 
Table Ixl However, only two of these operators are new and nontrivial. 

It is straightforward to show that the other three spurion combinations in (|A26|) do not generate any additional 
0(a 4 ) operators, so we choose to omit their discussion here. All that remains is the somewhat tedious task of replacing 
the taste spurions with specific taste matrices. There are three non-trivial pairings of tastes, leading to three different 
taste matrix substitutions: 



[V, Ax P] with [V, Ax P] => F = F' = £ 5 

[V, Ax P] with [V, A x T] => F = F' = ^ 

[V, AxT] with [V, A x T] => F = F' = £ pa 



(A27) 
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Operator 


Keep f 


69. 


Str(F L T,F R T 1 1 F^EF^E') Yes - 


combined with 74 


70. 


Str(Ft EFr FrE 1 F L ) 


No 


71. 


Str(F L F£SF fl FBE t ) 


No 


72. 


StT(F L FlT,F R F R ^) 


No 


73. 


Str(F L T,F' R F R E i F L ) 


No 


74. 


Str (F L EF' R ^F' L EF R ^) Yes - 


combined with 69 


75. 


Str(F L E J F fl E t )Str(F^EF^E t ) 


Yes 


76. 


Str(FiFl)Str(F fl Fi) 


No 


77. 


Str(F L EFi ? E t )Str(Fi,EF Ji E t ) 


Yes 



TABLE IX: The nine linearly-independent 0(a 4 ) operators composed of taste spurions Fl, F r , F' l and F' B 



Operator Keep? 

78. Str{F L F L HF l R F' R p)+p.c. No 

79. Str(F L EF^E t F L EF^E t ) + p.c. Yes 

80. Str(F L F L )Str(F' R F' R )+p.c. No 

81. Str(Fi,EF^E t )Str(F L EF^E t ) +p.c. Yes 



TABLE X: The four linearly-independent 0(a i ) operators composed of two Fl and two F R taste spurions. 



Note that, in the case of taste T with T, the two pairs of taste indices must be contracted separately because S 6 
respects both S'0(4)-taste and Lorentz symmetry. The three corresponding sets of chiral operators are given in 
Table IXXl 



b. Spin S, P, or T four-fermion operators with spin S, P, or T four-fermion operators 

The taste spurions corresponding to operators with spins S, P, and T are Fl and Fr, and can be combined in five 
parity-invariant ways: 

(FlFr^Fl'Fr) 

{V l ¥ l )Cf l 'f'r)+p.c. 

(F l Fr)(F l 'K')+P.c. 

(FlFl)(Fr'Fr')+ P .c.. (A28) 

In this case we build operators out of the left-handed basis defined in i|A12(l . Neither FrJ^ nor EFr is supertraceless, 
so each transforms as a bifundamental under the chiral symmetry group. Thus their product contains slightly more 
singlet irreps, and therefore generates more operators, than in the spin V case. 

Operators of type (FlFr)(Fl Fr ) correspond to singlet irreps in the product of four bifundamentals, i.e. the 
product of Fig. 13 with itself. There are twenty-four singlet irreps in this product, but we only show, in Table [XTI the 
three that lead to new 0(a 4 ) operators. 

In operators of type ((FlFl)(Fl Fr ) + p.c.) and ({FlFr)(Fl Fl ) + p.c), one pair of taste spurions is identical. 
Thus the number of such operators equals the number of singlet irreps in the product of one symmetric product of 
two bifundamentals with another two bifundamentals, i.e. the product of Fig. 01 with Fig. This product contains 
fourteen singlets, so there are twenty-eight such operators. However, only eight are new, and these are listed in 
Table IXill 

Finally, in operators of type ((FlFl)(Fl Fl ) + p.c.) and ((FlFl)(Fr Fr ) +p.c), both pairs of taste spurions 
are identical. Therefore these operators correspond to chiral singlets in the product of one symmetric product of 
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Operator Keep? 

82. Str(fiE t Fl'E t )Str(Ej5s^') Yes 

83. Str(JiE t ^'E t )Str(E?5)Str(EF^') +p.c. Yes 

84. Str(F t E t )Str(F|E t )Str(EF Ji )Str(EF^) Yes 



TABLE XI: The three non-trivial 0(a 4 ) operators composed of the four taste spurions Ft,, Fr, F' l and F' R 





Operator 


Keep? 


85. 


Str^E^E 1 " Fi,E t )Str(E^) + p.c. 


Yes 


86. 


Str^Et^Et F t E t )Str(EF fl ) +p.c. 


Yes 


87. 


Str(F i E t F i E t )Str(Fi,E t )Str(E^) + p.c. 


Yes 


88. 


Str(F^E t Fi,E t )Str(F i E t )Str(E? ; fl) + p.c. 


Yes 


89. 


Str(F L E t Fi,E t )Str(_FiE t )Str(EF^) + p.c. 


Yes 


90. 


Str(F^E t F L E t )Str(Fi,E t )Str(EF fl ) + p.c. 


Yes 


91. 


Str(F z ,E t )Str(Fi,E t )Str(Fi,E t )Str(EF^) + p.c. 


Yes 


92. 


Str(F{,E t )Str(Fl,E t )Str(FiE t )Str(EF R ) + p.c. 


Yes 



TABLE XII: The eight non-trivial 0(a 4 ) operators composed of two Fl, one F' L and one F' R taste spurions. 

bifundamentals with another symmetric product of bifundamentals, or of Fig. 0] with itself. There are ten such 
singlets, and therefore twenty such operators, but only fourteen are new. We list these in Table IXlTTl 

As we know from the previous section, spin T four-fermion operators generate a subset of those chiral operators 
generated by spins S and P; we can therefore ignore them. Thus, when pairing four-fermion operators, there are 
three combinations requiring different taste matrix substitutions: 

[S, PxV] with [S, P x V] =*> F = ^, F' = 
[S, PxV] with [S, P x A] F = F' = U 

[S, PxA] with [S, P x A] F = fa, F' = U ■ (A29) 
In all of these operators /i and v must be separately contracted to maintain 5 l O(4)-taste and Lorentz symmetry. 





Operator 


Keep? 


93. 


Str(htf Fl^F^F^) +p.c. 


Yes 


94. 


Str(Fx,E^E^ iE tfv E t) +p .c. 


Yes 


95. 


StT(F L ^F L ^F^)Str(F^)+p.c. 


Yes 


96. 


StT(F' L ^F' L ^F L ^)Str(F L ^)+p.c. 


Yes 


97. 


Str(F L E t F I ,E t )Str(F|E t F],E t ) + p.c. 


Yes 


98. 


Str^E+Fi.E^StrOFiEt.FlEt) + p. c . 


Yes 


99. 


Str^E^iE^Str^E^Str^Et) + p.c. 


Yes 


100. 


Str(F Zj E t Fi,E t )Str(F i E t )Str(JXE t ) + p.c. 


Yes 


101. 


Str(^XE t ^X2 t )Str(?iE t )Str(FLE t ) + p.c. 


Yes 


102. 


Str(F L E t )Str(Fz,E t )Str(FlE t )Str(Fl,E t ) + p.c. 


Yes 


103. 


Str(F L E t F I ,E t )Str(EF^E^) + p.c. 


Yes 


104. 


Str(F L E t Fz,E t )Str(EF^)Str(EF^) +p.c. 


Yes 


105. 


Str(F L E t )Str(F I ,E t )Str(EF]jEF^) +p.c. 


Yes 


106. 


Str(F L E t )Str(F L E t )Str(EF^)Str(EF^) + p.c. 


Yes 



TABLE XIII: The fourteen non-trivial 0(a 4 ) operators composed of two Fl and either two F' L or two F' R taste spurions. 
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Operator Keep ? 

107. Str(F L T,F R T,F R T, t F L T, 1: ) Yel 

108. Str(F i E^EfflE t )Str(fiE t ) Yes - combined with 109 

109. Str(F L E J F R E t FlE t )Str(EF^) Yes - combined with 108 

110. Str(F L E J F R E 1 ')Str(fiE t )Str(EF^) Yes 

111. Str(fiE t FlE t F L EF fl E t ) + p.c. Yes 

112. Str(^E t ^E t )Str(F t EF R E 1 ') + p.c. Yes 

113. Str^EtFLEFjjE^Str^E 1 ") + p.c. Yes 

114. Str(^E t )Str(^E t )Str(_FiEF fl E t ) +p.c. Yes 



TABLE XIV: The eight non-trivial 0(a 4 ) operators composed of spurions Fl, Fr and two F's . 
Table IXXll contains all of the corresponding operators in the chiral Lagrangian. 

c. Spin V or A four-fermion operators with spin S, P, or T four-fermion operators 

Operators that arise from the insertion of one operator with spin V or A and one with spin S, P, or T necessarily 
have mixed tastes. Each four-fermion operator contributes two taste spurions, so the corresponding chiral operators 
contain two -F's and two F's. The five parity-invariant ways of combining these spurions are 

(F l F r )(F l Fr) 
{F l F l )(FlF r ) + p.c. 
{F L F R )(F L F L )+p.c. 
{F l F l ){Fl¥ l )+p.c. 

(F L F L )(F R F R )+p.c. (A30) 

We do not need to use primes, since the spurions from the two four-fermion operators are already distinguished by the 
presence or absence of the tilde. As before, we build our operators using a left-handed basis, and in this case we must 
use both taste spurion bases defined in i|A7|l and (|A12(I . Recall that under chiral symmetry transformations, Fl and 
S-FrS^ transform as adjoints, while FlY$ and £Fr transform as bifundamentals. Thus the group theory counting 
will be slightly different than in the two previous sections, and each of the five spurion combinations in (|A30|) must 
be considered separately. 

First consider operators of type {Fl,F r )(FlF r ) . These come from the singlet irreps in the product of two adjoints 
with two bifundamentals. There are fourteen, but only four correspond to new NLO operators. Next consider the 
third type, (FlF r )(F[ j F[ j ) + p.c Because there are two F^'s, chiral operators correspond to chiral singlets in the 
product of two adjoints with the symmetric product of two bifundamentals. This product contains eight singlets, but 
only leads to four new NLO operators. Finally, the other three combinations in (|A30(I can be shown not to produce 
any new operators, so we do not discuss them further. Table IXIVI shows the eight new operators that can be built 
out of two F's and two F's. 

In this mixed-taste and mixed-spurion case, there are four possible four-fermion operator combinations and corre- 
sponding taste matrix substitutions: 

[V, A xP] with [S,PxV] => F = £ 5 , F = ^ 

[V, Ax P] with [S,Px A] => F = £ 5 , F = ^ 5 

[V,Ax T] with [S,Px V] => F = ^ u , F = f p 

[V, AxT] with [S, P x A] => F = ^u, F = £ p5 . (A31) 

Once again, all index pairs are contracted separately, thereby maintaining S'0(4)-taste and Lorentz symmetry. The 
chiral operators of C(a 4 ) arising from the insertion of a spin V or A four-fermion operator and the insertion of a spin 
S, P, or T four-fermion operator are given in Table IXXIII 
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d. Double insertions involving four-fermion operators in S 6 

In this section we determine mesonic operators of C(a 4 ) which arise from the double insertion of four-fermion 
operators, at least one of which is from S 6 . As in the single four-fermion operator insertion case, the chiral 

FF(A) 

structures are identical to those for insertions of the corresponding operators in S 6 , but new index structures 
arise which break the 50(4) taste symmetry. Also like before, we are not interested in operators which do not 
break this symmetry, since these have already been completely determined. Thus we first consider insertions of two 

F F( B) 

four-fermion operators from S 6 . After this discussion, it will be easy to see that combinations of one FF(A) 
operator with one FF(B) operator cannot generate anything additional, and in fact only generate operators that 
respect SO(4)-taste and Lorentz symmetries. 

Although the spin order is reversed as compared to previous sections, we begin with double insertions of [T„ x V M , A^]. 
This is because it is easy to guess the correct SO(4)-taste symmetry breaking chiral operators. They have the same 
chiral structure as a subset of operators that come from double insertions of [£, PxV,A], namely those which arise from 
L — R cross-terms (operators 82-84). A quick glance at Table IXXll reveals that there is only one possible symmetry- 
breaking index contraction: p = v. As this contraction does not break any of the discrete lattice symmetries, it must 
be present at some order in our power counting. We can show, however, that it appears at 0(a 4 ) using an argument 
similar to those in Section [A 31 For fixed p and //, the general spin-taste index structure of these operators (ignoring 
factors of £5 in A^) is the following: 

^ w ® aOC^A* ® ® m'X^V ® aO (A32) 

Because mesonic operators of 0(a 4 ) have no derivatives, invariance under the lattice symmetries implies that they 
must be singlets under hypercubic rotations. As was discussed previously, both pp and p'p' contain a singlet and a 
3-d irrep. Thus their product contains two overall singlets - one from the product of the singlets and one from the 
product of the 3-d irreps. Only the latter gives new operators, for it contains a part in which p = p' , confirming our 
initial guess. The £0(4)-taste symmetry breaking operators arising from double insertions of [T M x V^, A^] are given 
in the first two panels of Table IXXIIII 

Next we proceed to mixed-spin operators, i.e. the result of an insertion of [V^,A^ x T M ] with an insertion of 
[Tf, x Vfj,, Afj]. These operators have identical chiral structures to the ones in the top portion of Table IXXllI 13 Once 
again, there is only one possible £0(4)-taste breaking index contraction, p = p, but we must still show that this is the 
correct mapping of the four-fermion operator product. 14 In this case, the operators have the following index structure 
at fixed p and p': 

X! ® A^X/-* ® vfJ-)(n'v' ® V ® A 1 ') • (A33) 

Just as before, one hypercubic singlet in this product comes from the two 3-d irreps, and contains a part with p = p' . 
This generates the 0(a 4 ) mixed-spin operators in the third panel of Table IXXIIII 

Finally we address double insertions of [V^ , A^ x T M ] . This will produce operators with the same chiral structure 
as double insertions of [V, AxT]. Such operators are given in the bottom portion of Table IXXl In this case, multiple 
combinations of index contractions are possible: fi = p and p = p, v = a. Both are allowed by the lattice symmetries, 
but we must determine whether either arises at C(a 4 ). After fixing p and p! ', the index structure is 

^ ^ (M® A^Xm ® vn){jx' ® p'v')(p' ® v'p!) . (A34) 

Now we can see that the only nontrivial way of producing a hypercubic singlet is p = p' , just like in the two previous 
examples. The indices v and v' are independently summed from the beginning. We conclude that only the first of 
the two possible index contractions arises from the double insertion of four fermion operators. The final panel of 
Table IXXIIII shows the 0(a 4 ) chiral operators generated by double insertions of [V^,A M x T)J. 

After these examples, it is easy to see that the combination of an FF(A) four-fermion operator with an FF(B) 
four-fermion operator cannot generate any additional operators. This is because, in the ^^(^4) operator, spin and 



Note that we do not need the " — fiu5 " parts of [V M , xTJ, since they can be removed by adding an appropriate amount of [V, AxT], 

as discussed in Sec. IA 31 This simplifies the final taste matrix substitution. 

Note that v = p does not produce different operators since these are all just dummy indices. 
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taste indices are completely uncorrelated. Thus these operators contain only a singlet under hypercubic rotations, 
and can only form a overall singlet by combining with the singlet component of the FF(B) operator. But doing so 
removes the correlations between taste indices. Such mixed-symmetry double insertions only modify the coefficients 
of the 0(a 4 ) FF(A) operators. 

A natural question that arises is what happened to the second possible set of index contractions in the \V^, x T^] 2 
case, in which both /i = // and v = v' . This would lead to mesonic operators such as 

{Str^E^E+^E^Et) +p.c.} , (A35) 

i.e. with both /i and v appearing four times. Such operators are consistent with all of the lattice symmetries, so they 
must be present in the chiral Lagrangian at some order. The answer is that these operators are present, but they 
are of 0(a 8 ). It is clear from the previous examples that contracting the spin indices to make hypercube singlets can 
never contract all of the taste indices, some of which are independently summed. Thus chiral operators such as (|A35|) 
can only arise if all of the tensor indices are coupled in the quark-level operator itself. This requires a local eight-quark 
operator, e.g. 

X) 2(5(7/1 ® $M\Q^ ® U)Q][5(7„ ® C)Q][Q(7m ® U)Q] ■ (A36) 
which enters only at 0(a s ). 

A more physical way of looking at this is in terms of gluon exchange. Taste-changing interactions are due to 
the exchange of hard gluons with momentum k ~ ir/a. To obtain a complete correlation of taste indices requires a 
complete correlation of the hard gluon momenta. This can arise directly, for example, with four bilinears connected by 
four gluons joining at a four-gluon vertex, as in the above operator. The presence of four hard gluon propagators gives 
the factor (l/k 2 ) 4 ~ a 8 . In contrast, the 0(a 4 ) mesonic operators that we are considering arise from two separate 
insertions of 0(a 2 ) four-fermion operators tied together by gluons from the O(a ) Lagrangian. 15 Thus these gluons 
do not violate taste symmetry. The correlation between taste indices in the two four-fermion operators can only come 
indirectly, through the Lorentz structure. 

We can check the conclusion that operators such as (|A35|) do not appear at C(a 4 ) by examining the counterterms 
generated by other S'0(4)-taste symmetry breaking operators. Consider four-pion vertices generated using the 0(a 2 p 2 ) 
operators from the top portion of Table IXIXl Form a single diagram with only four external pions by joining the four 
pion legs with derivatives to form a loop. The loop integral couples the indices on the derivatives and gives rise to a 
quartically divergent counterterm with the following four tastes of external pions: fj,v, i/fj,, fia, and er/i. This is of the 
same form as resulting from the 0(a 4 ) operator that we derived above. 

The general conclusion we draw from this discussion is that one cannot determine the staggered chiral Lagrangian 
simply by testing the invariance of all possible operators under all of the lattice symmetries. While this procedure 
does filter out incorrect operators, it does not lead to the correct power- counting. The only foolproof method is 
therefore to map each quark-level operator directly onto chiral operators. 

5. Additional source terms for vector and axial currents 

Up to this stage we have not explicitly included source terms for left- and right-handed (or equivalently vector and 
axial) currents. The standard method for doing so is to promote derivatives to covariant derivatives, 

— » £>„E = M E - tf M E + »Er„ , — » D^Y? = - + »E% , (A37) 

and enforce local chiral symmetry |42l |43| . Here and r^ are, respectively, matrix sources for the left and right- 
handed currents. One also needs the field strengths associated with these sources, from which one can build contact 
terms at NLO. These do not contribute to S-matrix elements or decay constants, however, so we do not list them 
here. 



15 These lead to a local operator in the chiral Lagrangian because the gluons are confined to a small region by confinement. In effect, this 
replaces four powers of a in the eight fermion operator with i? 4 , where R ~ 1/Aqcjj is the confinement length scale. 
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Operator 


Keep f 


115. 


StWD FtFtY.D ^"i i „ r 
o bi \i-J l z-ilj^ z_* j -\- p.t^. 


Yes 


116. 


StT(D„F L T,D^ F L ) + p.c. 


Sometimes 


117. 


Str(D M F i EF flJ D M E t ) +p.c. 


Yes 


118. 


StT(D„F L ED^T,F R ^) +p.c. 


Sometimes 



TABLE XV: The four a 2 source operators with one D^F and one Lie derivative. Note that the covariant derivative only acts 
on the object directly following. The operators which are to be kept "sometimes" can be removed using the equations of motion 
if the underlying operator is of type FF(A). 

In continuum %PT, the NLO Lagrangian takes the same form as that without the sources except for the change 
of regular to covariant derivatives, aside from aforementioned contact terms. Possible terms involving covariant 
derivatives acting on M or M\ for example, can be removed using the equations of motion (EoM). This turns out 
not to be true for S^PT. There are additional terms involving covariant derivatives acting on the taste spurions which 
cannot be removed, and which play an important role in the breaking of 5*0(4) for axial decay constants. Of course, 
in addition to these terms one must change regular to covariant derivatives in all the 0(a 2 p 2 ) terms given previously. 

The extra "source terms" are built from a covariant derivative acting on a taste spurion, a second taste spurion 
(since they always come in pairs), and a Lie derivative (to match the Lorentz indices). Thus they are of 0(a 2 p 2 ) and 
are produced by single insertions of four-fcrmion operators. We run through the possibilities in turn. 

a. Single insertion of FF(A) operators 

We will go through the group theory and operator enumeration rather quickly, as the process is the same as in 
all previous sections. We will concentrate instead on eliminating operators through use of the EoM, as this is where 
there is a difference between staggered and continuum %PT. 

Operators arising from a single insertion of [V, A x F] can be built out of the left-handed basis in l|A7|) along with 
two additional elements: 

D^F L , ZD^Fr^ . (A38) 

Because Fl and Fr transform differently than E and T,\ the covariant derivative acts in the following way on the 
taste spurions: 

D^F L = O^Fl - i %, F L ] , D^F R = d li F It -i [r M , F R ] . ( A39) 

Once Fl^r are set to fixed taste matrices, the partial derivative pieces will vanish. However, the commutator terms 
can still contribute to the left- and right-handed currents, and therefore to PGB decay constants. 

It is clear from (|A39(I that taste S cannot contribute, since then the commutator vanishes. Thus we need only 
consider tastes P and T, in which case Fl and Fr are supertraceless, and transform as adjoints under the chiral 
symmetry group. Moreover, D^Fl and D^Fr are also supertraceless because of the cyclicity of the supertrace, so 
they too are adjoints. Thus operators with one taste spurion, one covariant derivative acting on a taste spurion, and 
one Lie derivative will correspond to the singlet irreps in the product of three adjoints. There are two singlets in this 
product, and therefore four operators. These operators are given in Table IXVl 

We can remove some of these operators using the EoM. First recall (qualitatively) how this is done in the case 
of covariant derivatives acting on mass spurions. One such operator is Str(D^MD^T,^). Integrating by parts, the 
covariant derivative can be moved onto the field, resulting in — Str(MD fl D fl 'E^). Then D^D^Y^ can be replaced with 
operators involving either one or no derivatives acting on each of E or E^ using the EoM, all of which are already 
included in the standard NLO effective Lagrangian. By trying to follow the same procedure for operators involving 
taste spurions, we will see where it breaks down. 

For example, first consider operator 115. We must be careful when integrating by parts because there are more 
than just two terms in this operator - we cannot simply move from one piece to the other and add a minus 
sign. Derivatives can be moved around using integration by parts because we are integrating over all spacetime in the 
action, so total derivative terms become surface terms and can be neglected. We therefore consider a total derivative 
term which gives rise to operator 115: 



D^FlFlHD^) = (D^F L )F L ZD^ + F L (D^F L )HD^ + F L F L D^D^ + F L F L X(D „D M £t) . (A40) 
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Operator 


Keep f 


119. 


Str(D ll F L T, J F L D ll ^ J ) + p.c. 


Yes 


120. 


Str(_D M fz,-D M £ T -FL£ T ) +p.c. 


Sometimes 


121. 


Str(D M Fz,E T )Str(Fz,D M E T ) +p.c. 


Yes 


122. 


Str^D^FiD^STjStr^LE 1 ) +p.c. 


Sometimes 


123. 


StrfD Fr PuYD Y^) 4- n 
o\>v\±J n Z-iU ^Z-i j -t p.c. 


I CO 


124. 


Sti(D^¥ L D^Y,F R ) + p.c. 


No 


125. 


Str( J D M JYE t )Str(EF^ED M E t ) + p.c. 


Yes 


126. 


Str( J D M ?i J D M E t )Str(EF^) +p.c. 


Sometimes 



TABLE XVI: The eight a 2 source operators with one D^F and one Lie derivative. Notation as in Table IXVl 

After integrating over all spacetime and taking the supertrace, the l.h.s. vanishes, so the four operators on the r.h.s. 
are all related. The first is operator 115, while the second is operator 116. The third operator is from the LO 
Lagrangian, while the fourth can be related to other standard operators using the EoM. We learn two lessons from 
this exercise: first, that operators 115 and 116 are redundant - they are related by integration by parts; second, that 
we cannot remove both using the EoM — and we choose to keep operator 115. The difference from the example given 
above of the covariant derivative acting on the mass matrix is the presence of two taste spurions in each operator. 
Integration by parts can never remove the covariant derivative from both taste matrices and produce something such 
as D^D^T,, which can be removed using EoM. The same holds for operators 117 and 118; we can remove only one of 
these two using the EoM, and we choose to keep operator 117. 

Operators arising from a single insertion of [S, PxF] can be built similarly out of the left-handed basis in (|A12|) 
plus the elements 

D^FlYJ, VD^Fr. (A41) 

In this case Fl and Fr transform in the same way as S and Yt\ so 

D^Fl = d^F L - U^Fl + iF L r^ , D^Fr = B^Fr - ir^FR + iFr£^ . (A42) 

Neither F^Y^ nor D^FlYJ (or their parity conjugates) are supertraceless, so they transform as bifundamentals. Thus 
all operators with one F, one D^F, and a Lie derivative come from singlet irreps in the product of two bifundamentals 
with one adjoint. There are four such singlets, and eight such operators, given in Table IXVll 

We can now eliminate operators using integration by parts and EoM. Operator pairs 119 & 120, 121 & 122, 123 & 
124, and 124 & 125 are each related in this way, and we choose to keep only the first operator in each pair. In fact, 
operator 124 becomes equal to operator 123 once one uses the result that D^FlFr) — 0, i.e. that the spurions are 
momentum-independent and their product has no taste. Thus we are left with four additional operators. 

As in the previous sections, four fermion operators with spin T do not lead to any additional mesonic operators, 
and simply modify some of the unknown coefficients of the spin S and P operators. 

Inserting the appropriate taste matrices is now straightforward. We simplify expressions using the fact that d^F = 0. 
We also write operators in terms of £^ and r M because the covariant derivative acts differently on the two types of 
taste spurions. The final "source term" operators are listed in Table IXXIVI 

b. Single insertion of FF(B) operators 

In Section |A 31 we discussed in detail how to map onto 0(a 2 p 2 ) mesonic operators arising from a single insertion 
of a four fermion operator in 5 6 , so we just summarize the result here: the chiral structures are the same as 
those of single FF(A) insertions, but the index structures are different, breaking 5'0(4)-taste and Lorentz symmetry. 
Simply put, each FF(B) operator is identical to the corresponding FF(A) operator, but with the derivative indices 
contracted with a pair of taste indices. However, there is one important difference: because the indices on the covariant 
derivative and the taste matrices are correlated, the number of operators can no longer be reduced using the EoM. 
Thus there are more a 2 source operators corresponding to a single insertion of [V^ x T^] or [A^ x T^] than those those 
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from a single insertion of [V, A x T]. 16 There are, however, the same number of a 2 source operators from [T M x V^, A^] 
as from [S x V, A] because the doubling from the inability to use the EoM is counteracted by a halving because there 
are only L — R cross-terms. Table IXXVI lists all of the resulting FF{B) "source term" operators. 

Note that [V^ x T M ] only generates single supertrace operators, so there are no "source term" hairpins with tensor 
taste. This has interesting consequences for 50(4)-taste symmetry breaking in PGB decay constants. 
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Generic 


Specific 0(a 


2 p 2 ) Operator 


Operator 


[V, Ax P] 


[V,AxT] 


1 + 3 


StrO^Et^EEt&Ses) +p.c 


StrO^Et^EE^pE^) + p . c . 


2 




Str(a M Et^ p a M Ef p „) 


6 


Str(9 M E t a M E)Str(C 5 E^ 5 E t ) 


Str(9 M E t 9 M E)Str(^ P E^E t ) 


7 


Str(Ea M E t ^ 5 )Str(E t a M EC5) 


Str(Ea M E t ^p)Str(E t 9 M Ee p ,) 


10 


Str(E9 M E t e5E9 M E t e 5 )+p.c. 


Str(Sa M E t ^ p E,9 fl E t Cp,) +p.c. 


13 


Str(Ea M E t C 5 )Str(E9 A1 E t C 5 ) + P-c. 


Str(Ea Al E t ^ p )Str(E9 p E t Cp.) + P-c 



Generic 


Specific 0(a 


2 m) Operator 


Operator 


[V, AxP] 


[V,AxT] 


15 + 22 


Str(fcE&Af + )+p.& 


Str^E^pAft) +p.c. 


16 + 21 


StrfeEAftEesE 1 ") + p.c. 


Str(e M „EM t E^pE t ) + p.c. 


17 + 23 


Str(C 5 Ee5E t )Str(EM t ) +p.c. 


Str(^„E^ Al E t )Str(EM t ) +p.c. 



TABLE XVII: Operators in the staggered chiral Lagrangian arising from a single insertion of an S 6 operator with spin V 
or A. Implicit sums over repeated indices follow the summation convention defined in Q, and p.c. indicates parity conjugate. 
Each chiral operator has an independent undetermined coefficient of 0(a 2 ). Partial derivatives act by convention only on the 
object immediately to their right. To include left- and right-handed sources, derivatives should be replaced with covariant 
derivatives. M and serve as (pseudo) scalar sources; if such sources are absent then M = — M. 



Note that the operators resulting from generic forms 115 and 116 are the same up to a sign, so the doubling of operators does not apply 
in this case. 
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Generic 


Specific 


0(a 2 p 2 ) Operator 


Operator 


[S, PxV] 


[S, Px A] 


36 




Str(E9 M E T ^ 5 E T a M E$ 5 „) 


38 + 39 


Str(9 M E T 9 M EE T ^)Str(E^) +p.c. 


Str(a M E T 9 M EE T ^ 5 )Str(E^ 5y ) +p.c. 


41 


Str(<9 M E T £„ ) Str (d,, E£„ ) 


Str(a M E T ^5)Str(a M E^) 


42 


Str^a^E'^EjStr^^E 1 )Str(E£„) 


Str^dpE'S^EjStr^sE 1 )Str(E^ 5 ^) 


43 


Str(9 M E9 M E T ^E T ^E T ) + p.c. 


StrtcVScVE^sE^E 1 ) +p.c. 


44 


Str(d M E'£,,cVE '<<;,,) + p.c. 


Str(5 M E T ^ 5 a M E T ^ 5l ,) + p.c. 


45 


^tr(r) Vfrl Y.^trtf V,V 5Th 4- n r 


Strfr) T^r) yi'lStrf/ 7 =V,Vc Vh + w c 

Obl^C^^j O LI ^5^^ ^ ^ p.C 


46 


Str(a M Ed M E t C,E t )Str(C„E t ) + p.c 


Str^EcVE^sE^Str^E 1 ") + p.c. 


47 


Str^E^) 2 + p.c. 


Str(a M E t e„5)Str(9 fl E t C 5 .) + p.c. 


48 


Str(9 M E t 5 M E)Str(^E t ) 2 + p.c. 


Str(9 M E t 9 M E)Str(C,5E t )Str(C5,E t ) + p.c. 



Generic 


Specific 0(a 


2 to) Operator 


Operator 


[S,PxV] 


[S, P x A] 


52 + 59 


Str(C M M t )Str(E^ M ) 


Str(C M5 M t )Str(EC 5 M) 


53 + 58 


Str(C M E t )Str(EC M EAf t ) +p.c. 


Str(^ 5 E t )Str(EC 5 MSA/ t ) +p.c. 


54 + 60 


Str(^E t )Str(Ee M )Str(EM t ) + p.c. 


Str(e M 5E t )Str(Ee 5M )Str(EM t ) + p.c. 


61 


Str(^Et^Mt) +p.c. 


Str(^ 5 Et^Mt) + p.c 


62 


Str^E^E^StrtEM 1 ") +p.c. 


Str(^ 5 E t C 5 MS t )Str(EM t ) +p.c. 


63 


Str(^ M E t )Str(e M A/ t ) + p.c. 


Str^sE^Str^Mt) +p.c. 


64 


Str(£ fl E t )Str(£ M E t )Str(EM t )+ p.c. 


Str(^5E t )Str(^ 5M E t )Str(EM t ) +p.c. 


65 


Str^Et^EtME 1 -) +p.c. 


Str(e Al5 Et^E t MEt)+p.c. 


66 


Str(^E t ^E t )Str(ME t ) +p.c. 


Str(e M5 E+^ M Et)Str(MEt) +p.c. 


67 


Str^E^Str^EtME 1 ") + p.c. 


Str^sE^Str^StME 1 ") +p.c. 


68 


Str(^E t )Str(^E t )Str(ME t ) + p.c. 


Str(^ 5 E t )Str(e 5 MS t )Str(ME t ) +p.c. 



TABLE XVIII: Operators in the staggered chiral Lagrangian arising from a single insertion of an S 6 operator with spin S 
or P. Notation as in table lX*VIll 
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Generic 


Specific 0(a 2 p 2 ) Operator 


Operator 


[V M x T M ] and [A^ x T^] 


1 + 3 


Str(9 M E T a M EE T f MI/ E£„ M ) +p.c. 


2 


Str(a M E T 5 M „a M E^ M ) 


5 




7 


Str(E9 M E t ^,)Str(E t 9 M Ee^) 


10 


Str(Ea M Ete M ,E9 M Et^ M )+p.c. 


13 


Str(Ea M E t ^,)Str(Ea M E t C„ M ) + p.c. 



Generic 


Specific 0(a 


2 p 2 ) Operator 


Operator 
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Str(E9 M Et^Eta M E^) 


Str(Ea M Et^ 5 Eta M Ee 5M ) 


38 + 39 


Str(a fl E t a fl EE t ^)Str(E^) +p.c. 


Str(a M E t 9 M EE t ^5)Str(EC5 M ) 


41 


Str(a M Et^ M )Str(a M E^) 


Str(a M Et^5)Str(a M EC 5M ) 


42 


Str(9 M E t a M E)Str(e M E t )Str(E^) 


Str(a M E t a M E)Str(C fl 5E t )Str(E^) 



TABLE XIX: Rotationally non-invariant operators in the staggered chiral Lagrangian arising from a single insertion of S t 
operators. There is an implicit summation over both fj, and v with the constraint v ^ fj,. 



Generic Op. 


[V, Ax P] with [V, Ax P] 


69 + 74 
75 


StrfeE^EtesECsE 1 ") 
StrfeECsE^StrfeEesE 1 ") 



Generic Op. 


[V, Ax P] with [V, A x T] 


69 + 74 


StrC^ECsE^E^E 1 ') +pc 


75 


StrfeE^E^Str^E^E 1 ) 


77 


Str(f 5 E^Et)Str(^E^Et) 


79 


Str^E^EtCsE^E 1 ") +p.c. 


81 


Str(f 5 E^„E t )Str(f 5 E^E t ) + p.c. 



Generic Op. 


[1/, ylxT] with [V, A x T] 


69 + 74 


Str^E^E^paE^pEt) +p.c. 


75 


Str(e M ,E^E t )Str(^ CT E^ p E t ) 


77 


Str^E^E^Str^pE^E 1 ") 


79 


Str(^„E^p CT E t ^ I/M E^ CT pE t ) +p.c. 


81 


Str(^ I/ E^ CT E t )Str(f„^Ef CTp E t )+p.c. 



TABLE XX: C(a 4 ) operators in the staggered chiral Lagrangian arising from two insertions of an S , ( f F< ' A ' operator with spin 
V or A. Notation as in table IX Vlfl 
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Generic Op. 


[S, P x V] with [S, P x V] 


[5, P x 4] with [5, P x A] 


82 


Str(^ M E T ^E T )Str(E^ M E^) 


Str(^5E T ^5E T )Str(E^5 M E^ 5l/ ) 


83 


Str(^ M E t ^E t )Str(E^)Str(E^) + p.c. 


Str(^ 5 E t ^5E t )Str(E5 5fl )Str(E5 5t -) + p.c. 


84 


Str(£ M E T )Str(^£ T )Str(££ M )Str(E£,,) 


Str(^ 5 E T )Str(^5E T )Str(EC5, l )Str(E^5„) 


85 


Str(^ M ET^E T ^ET)Str(E^) +p.c. 


Str(^ 5 E T C5 Al S T ^5S t )Str(E^5^) +P-C. 


87 


Str(^ M E t ^ M E t )Str(^E t )Str(E^) + p.c. 


Str(^ M5 E t ^S t )Str(^5E t )Str(EC5 I ,) + p.c. 


89 


Str(^ /1 E t ^E T )Str(^E T )Str(E^) +p.c. 


Str(^5E t ^ 5 ET)Str(^5 fl E T )Str(EC5 I ,) + p.c. 


91 


Str(^E t )Str(^ /J E t )Str(^E t )Str(E^) + p.c. 


Str(f M 5E t )Str(^5 M E t )Str(^5E t )Str(EC5 I ,) +P-c. 


93 


Str(£ M E T £ M E T £„E T £„E T ) +p.c. 


Str(^ M5 E T ^ 5M E T ^ 5 E T ^5^E T ) +p.c. 


94 


Str(£ fl E t £„£ t £ fl E t £„£ t ) +p.c. 


Str(£ (j5 E t £ I ,5E t £5 M £ t £5 l ,E 1 ') + p.c. 


95 


Str(^E T ^ M E T ^E T )Str(^E T ) + p.c. 


Str(^ M5 E T ^5 M E T ^5E T )Str(^ 5l ,E T ) +p.c. 


97 


Str^E^E^Str^E^E 1 ") +p.c. 


Str^sE^^E^Str^sE^E 1 ") +p.c. 


98 


Str^E^E'OStr^E^E 1 ") + p . c . 


Str^sEt^sE^Str^sMS^s^E 1 ") +p.c. 


99 


Str(£ M E t £ fl E t )Str(£ I ,E t )Str(£ l ,E t ) + p.c. 


Str^sE^SMS^Str^gS^Str^s^E 1 ") + p.c. 


100 


Str^E^E^Str^E^Str^E 1 ") + p.c. 


Str^sE^sE^Str^E^Str^E 1 ") +p.c. 


102 


Strf£ Vflqt,.^ V'hqtrl'A y^QWt vtl i „ r 


Oll^^i5Zj ^OLI^S^Zj JOLI^iy5Zj ^011^5^^ j p.c. 


103 


Str(e M Ete M Et)Str(Ee,E^) 


Str(^ 5 Et^ M Et)Str(E^ 5 Ee 5 ,) 


104 


Str(C M E t ^S t )Str(E^)Str(E^) + p.c. 


Str(^5E t C5 M E t )Str(E^ 5 )Str(EC5.) + p.c. 


106 


Str(e M E t )Str(5 fl E t )Str(E^)Str(E^) 


Str(^5E t )Str(6 /J E t )Str(E^ 5 )Str(E^,) 



Generic Op. 


[S, P x V] with [S, P x A] 


82 


btr^E'^sE 1 )btr(E£ M E£ 5 „) 


83 


Strff„EV„«;E t )Str^Ef„ 1 )Str^Ef«i„ , ) + d c 


84 


Str(C M Et)Str(^ 5 Et)Str(E^)Str(EC 5 ,) 


85 


Str(C M E t ^E t ^5E t )Str(EC 5 ,) + P-c. 


86 


Str(^ 5 E t C5„E t ^E t )Str(EC, I ) + p.c. 


87 


Str(C M E t C M E t )Str(US t )Str(E^)+ p.c. 


88 


Str(e„sE t e 5! ,E t )Str(^E t )Str(E^)+p.c. 


89 


Str(C M E t C, 5 E t )Str(^E t )Str(Ee^) + p.c. 


90 


Str(C,5E t e M E t )Str(e 5 ,E t )Str(E^ 1 )+ p.c. 


91 


Str(^E t )Str(e M E t )Str(^ 5 E t )Str(E^)+ p.c. 


92 


Str(^ 5 E t )Str(C 5 „E t )Str(^E t )Str(E^)+ p.c. 


93 


Str^E^E^sSt^Et) +p . c . 


94 


Str^EtUS^MS^s.E 1 ") +p.c. 


95 


Str^Et^Et^sE^Str^E 1 ") + p.c. 


96 


Str (e,s E 1 "^. E 1 " ^ E 1 " ) Str(C M E 1 " )+ p.c. 


97 


Str^EtCME^Str^sS^s.E 1 ") + p.c. 


98 


Str(e M E t C,5E t )Str(C M E t e 5 „S t )+ p.c. 


99 


Sta(£ M Et£ M £ t )Sto(£„6E t )Str(£ 8 ,.Et) + p.c. 


100 


Str(^Et^ s Et)Str(^E+)Str(e Sl/ E t ) + p.c. 


101 


Str^sE^s.E 1 ") Str^E^Str^E 1 ") + p.c. 


102 


Str(^E t )Str(^E t )Str(^ 5 S t )Str(^E t )+ p.c. 


103 


Str(C M E t ^S t )Str(EUSC5.) + p.c. 


104 


Str(^E t C M E t )Str(E^5)Str(Ee 5 .)+ p.c. 


105 


Str(^ 5 E t e 5 ,E t )Str(EC M )Str(E^)+ p.c. 


106 


Str(^Et)Str(e M Et)Str(E^ 5 )Str(E^.) + p.c. 



TABLE XXI: 0(a 4 ) operators in the staggered chiral Lagrangian arising from two insertions of an Sq F ^ A ^ operator with spin 
S or P. Notation as in table IX VIll 
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Generic Op. 


[V, Ax P] with [S, P x V] 


[V, Ax P] with [5, P x A] 


107 


Str(£ 5 E£pE£ 5 E T £pE T ) 


Str(£ 5 E£p5E£ 5 E T £5pE T ) 


108 + 109 


Str(£ 5 E£pE£ 5 E t )Str(£pE t ) +p.c. 


StrfeECpsE^E^Str^spE 1 ') + p.c. 


110 


Str(£ 5 E£ 5 E T )Str(£pE T )Str(E£p) 


Str(^ 5 E^5E T )Str(^ 5 E T )Str(E^5A l ) 


111 


Strff EV E+fcE^cEh 4- n r 


^ L1 VsmS^ Sop^ So^So^ j p.c 


112 


Str^EtCME^Str^EesE 1 ") + p.c. 


Str(^5E t e5 M S t )Str(C 5 Ee5E t )+ p.c. 


113 


Str(e M Ete 5 EC 5 E t )Str(e M Et) +p.c. 


Str(^ 5 Ete 5 Ee 5 Et)Str(^Et) +p.c. 


114 


Str^E^Str^E^StrfeE^E 1 ") + p.c. 


Str^sE^Str^E^StrfeEesEt) + p.c. 



Generic Op. 


[V, AxT] with [S, P x V] 


[V, AxT] with [S, Px A] 


107 


Str^ECpE^E^Et) 


Str(ep.ECp 5 EC,pE t 6pE t ) 


108 + 109 


Str(C M ,EC P EC, M Et)Str(C P Et) +p.c. 


Str(^ECp 5 Ee, M Et)Str(e 5 pEt) + p.c. 


110 


Str E£„p E+ ) Str(^ p E t ) Str(E£ P ) 


Str(ep,E^pE t )Str(Cp5E t )Str(ECsp) 


111 


Str(^Et^Et^E^ M Et) +p.c. 


Str^psE^spE^p.S^pEt) + p.c. 


112 


Str^pE^pS^Str^E^pE 1 -) +p.c. 


Str(e pS Etf 5p Et)Str(^„E^ p Et) + p.c. 


113 


Str^pE^p.E^pE^Str^pE 1 ") +p.c. 


Str^psE^p.E^pE^StrfepEt) + p.c. 


114 


Str(f p Et)Str(^Et)Str(^ 1/ E^E+) + p.c. 


Str^psE^StrfepE^Str^E^pE 1 ") + p.c. 



TABLE XXII: 0(a 4 ) operators in the staggered chiral Lagrangian arising from two insertions of S FF operators, one with 
spin V or A, the other with spin S or P. Notation as in table IXVlTI 



Generic Op. 


[Tfi x V„] with [Tp x Vp] 


[Tp x Ap] with [Tp x Ap] 


82 


Str(^EtCpE t )Str(E^E?p) 


Str(^p5E t Cp 5 E t )Str(Ee5pSC5p) 


83 


Str(CpE^pE t )Str(ECp)Str(Eep) + p.c. 


Str(ep 5 E t ep 5 E t )Str(EC5p)Str(Ee 5 p) + p.c. 


84 


Str(C M E t )Str(CpE t )Str(ECp)Str(E^p) 


Str(Cp5E t )Str(^p 5 E t )Str(Ee5p)Str(EC5p) 



Generic Op. 


[Tp x Vp] with [Tp x Afj] 


82 


Str(CpE t Cp 5 S t )Str(ECpEC5p) 


83 


Str(CpE^p5E t )Str(ECp)Str(E^p) + p.c. 


84 


Str(epE t )Str(^p 5 E t )Str(ECp)Str(E6p) 



Generic Op. 


[Vp,^p x Tp] with [Tp x Vp] 


[V„,A„ x T M ] with [Tp x Ap] 


107 

108 + 109 
110 


Str(Cp,ECpE^pE t CpS t ) 
Str(^E^E^ M Et)Str(^Et) + p.c. 
Str(^,EC,pE t )Str(e M E t )Str(E^) 


Str^p.E^psE^pE^spE 1 ) 
Str(ep„Eep5E^pE t )Str(e 5 pE t )+p.c. 
Str(Cp,E^pE t )StrKp 5 E t )Str(EC 5 p) 



Generic Op. 


[Vp, Ap x Tp] with [Vp,Ap x Tp] 


69 + 74 


Str^p.E^pE^p^E^pE 1 ") + p.c. 


75 


Str(Cp,E^pE t )Str(Cp CT EC CT pE t ) 


77 


Str^p.ECp.E^Str^pE^pE 1 ) 


79 


Str(£p„E£p CT E t £„pE£ (T pE t ) +p.c. 


81 


StrCCp.E^E^Str^pE^pEt) + px . 



TABLE XXIII: C(a 4 ) operators in the staggered chiral Lagrangian arising from two insertions of Sg F operators. The indices 
fj, and v are separately summed, with the constraint that /i 7^ v. 
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Generic Op. 


[V, Ax P] 




[V, ixT] 


115 


iStr([£ M ,C5]e5SD M S t )+p.c. 


i Str([^ 


^]^SD,St)+ 2) .c. 


117 


tStr([4,,fc]E£ B £V,Et) +p.c. 


i Str([*„ 





Generic Op. 


[S, PxV] 


[5, PxA] 


119 


i Str([r M £„ - ^JE^Z^E) +p.c. 


i Str([r M &, B - ^ 5 £ M ]E£ Bl/ L> M E) +p.c. 


121 


i Str([r M £„ - ^^]E)Str(^Z3 (1 E) +p.c. 


i Str([r M ^5 - ^5^]E)Str(^-D M S) +p.c. 


123 


i Str([r M ^ -id^i^D^T,) +p.c. 


i Str([r M ^r, - (.^i^u^D^) +p.c. 


125 


i Str([r>£„ - ^]E)Str(£„.D M Et) +p.c. 


i Str([r M <e„ 5 - U^]X)Str{&„D^) + p.c. 



TABLE XXIV: 0{a 2 p 2 ) operators proportional to sources, arising from single insertions of S FF operators. 



Generic Op. 




115 


iSta^, +p.c. 


117 


i Str([^,£^]E£„ M .D M £ t ) +p.c. 


118 


i Str([^,^]S J D M S t E^ Al Et) + p.c. 



Generic Op. 


[T, x V,] 






123 


i Str([r^ ~ ^£^^ D^E) + p.c. 


i Str([r M £ M5 




125 


i Str([r M £ M - CM^]S)Str(C M D M E t ) + p.c. 


i Str([r>£ M5 - 


^ # ,]E)Str(£ B ,.A.E t )+p.c. 


126 


i Str([r M £ M ~ £ M £ M ]D M E)Str(E^) +P-C. 


i Str([r>£ M5 - 


C M 5^]^S)Str(E t e 5 M)+P-c. 



TABLE XXV: 0(a 2 p 2 ) operators proportional to sources, arising from single insertions of Sq F ^ operators. Indices are 
independently summed, with the constraint fi v. 



